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Chapter 1 

Foreword 



The Quantum Chromo-dynamics (QCD) is a theory of strong interactions 
— interactions between hadrons and, in particular, between their inner con- 
stituents. In QCD, the fundamental building blocks are quarks and gluons 
whose interactions are ultimately defined by the Lagrangian density 

>Cqcd = -1^;,^"''^'^ + {ij^'D^ - ruk) *^ (1.1) 

where denote the quark fields, 7''s the standard Dirac matrices, and 

= d^-igsf'Al (1.2) 
f;, ^ d^Al-d,Al + g,r'-AlAl, (1.3) 

where ^4^ are the gluon fields and Qs denotes the strong coupling constant. 
The matrices are the SU(3) generators and /"^'^ are the corresponding 
structure constants. The physics content of £qcd has turned out to be very 
rich, yet challenging to work out. The most rigorous approaches to probe the 
inner workings of QCD are the lattice simulations, which have demonstrated 
encouraging results e.g. for confinement and predicting the mass-hierarchy 
of light hadrons. The lattice-QCD, however, quickly meets its hmitations 
when the size of the studied system increases and it comes to describing scat- 
tering experiments. To apply QCD in such situation, perturbative methods 
to treat quarks and gluons are to be employed. The ultimate justification 
for the use of perturbative QCD (pQCD) tools lies in the fact that QCD 
enjoys what is known as asymptotic freedom — the strong interactions be- 
coming effectively weaker when the inherent momentum scale of the process 
is large, 2> 1 GeV^, or equivalently, when the probed distances are much 
smaller than the size of the hadron. As the strong interactions nevertheless 
bind the quarks and gluons, partons, together to make a hadron, the exact 
way they are distributed inside the hadrons cannot be neglected when ap- 
plying pQCD to hadronic colhsions. Intuitively, the structure of the hadron 
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should not, however, have anything to do with the collision, but is rather 
something that is inherent for the hadron itself. Prom the pQCD point of 
view, such property is known as factorization, and the relevant structure of 
the hadrons is encoded in parton distribution functions (PDFs) which are 
process-independent. In principle, the PDFs should be computable from 
^QCD but such task is far from being reahzed in practice any time soon. 
Instead, they must be inferred from various experiments with the help of 
pQCD — from global analyses. 

The role of the proton PFDs becomes emphasized in a hadron-hadron col- 
hder hke the CERN-LHC where the backrounds are often huge and the 
expected physics signals relatively weak. Interpreting the experimental mea- 
surements in a situation like this, requires reliable knowledge of the PDFs. 
Similarly, the detailed knowledge of the quark-gluon content of the bound 
nucleons is of vital importance in precision studies on the properties of 
the strongly interacting matter expected to be produced in ultrarelativistic 
Pb+Pb collisions at the LHC and e.g. Au+Au collisions at the BNL-RHIC. 

This thesis consists of two parts, the separate introductory part and the 
published four articles. The introduction begins by a technically detailed 
description of the DGLAP evolution — the pQCD-physics behind the global 
QCD analyses — as I understand it. I also discuss the fast numerical solving 
method for the DGLAP equations, which has been used in the numerical 
works of the published articles of this thesis. A write-up of the next-to- 
leading order (NLO) calculations for the deeply inelastic scattering (DIS) 
and the Drell-Yan (DY) dilepton production cross-sections, which are the 
data types that comprise most of the experimental input employed in the 
articles of this thesis, is also included. The formahsm of the inclusive sin- 
gle hadron production at NLO, the third type of experimental data utihzed 
in these articles, is described as well, although less rigorously. The intro- 
ductory part ends with a discussion of the global QCD analyses in general, 
with a special attention paid to the major work of this thesis [iVj, the NLO 
analysis of nuclear parton densities and their uncertainties. I have tried 
to avoid unnecessary overlap between the introductory part and the pub- 
hshed articles, but yet keep the introductory part such that it is logical and 
self-contained, without leaning too much on the published articles. The nec- 
essary background for understanding what is presented in this thesis is the 
basic knowledge of Quantum Field Theory and elementary phenomenology 
of High Energy Physics. 
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Chapter 2 

DGLAP evolution 



In this chapter, I will discuss the physics of parton evolution. Instead of only 
sketching general guidehnes, I will take a somewhat more detailed point of 
view, hoping this thesis would also serve as an elementary introduction to 
the subject. Much of what I present here can be learned from works of 
Dokshitzer et al. [lii2j and Altarelli [3j. 



2.1 Deeply inelastic scattering 




Figure 2.1: Schematic picture of the deeply inelastic scattering. The vari- 
ables k and P denote the incoming momenta, whereas k' and ki,. . . ,kn are 
all outgoing. The spin states are marked by s, s' and a. 

In a deeply inelastic scattering (DIS) a lepton projectile hits a target nu- 
cleon breaking it apart to a state X consisting of a plethora of various 
particles with invariant mass Mjr A/P, where M denotes the rest mass 
of the nucleon. In the simplest case the lepton is an electron or muon and 
the interaction is dominantly mediated by exchanging a virtual photon, as 
illustrated in Fig. O 
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In the target rest frame, the four-momenta of the particles can be chosen as 

(2.1) 



k = {E,k) ={E,0,0,E) 

k' = iE',k') = {E',E' sine cos (j),E' sine sin (l),E' cose) 

P = (P0,P) =(M, 0,0,0) 

q = {u,ci) =iE-E',k-k), 



where I have neglected the lepton mass. The standard invariant DIS- 
variables are 

Q2 = _g2 = 4EE' sin^ {e/ 2) 

X = ^-QL (22) 



y 



P ■ q V 

P-k ~ E' 



where the latter equahties refer to the target rest frame. The differential, 
spin-independent cross-section for this process can be written as 

where e is the electric couphng constant and 

L^"" = ^TV[/^'7^|^7'^] (2.4) 

(n,out|J^(0)|(P,cT),in)((P,a),in|jt(0)|n,out) 

are the leptonic and hadronic tensors. In contrast to the leptonic tensor L^^ , 
the non-perturbative nature of QCD makes it impossible to compute W^jy 
directly but its general form can nevertheless be written down without much 
further input. Indeed, since L^^ is symmetric under interchange of indices, 
the relevant part of the hadronic tensor should also satisfy W^y = Wy^,, 
and together with the definition ()2.5I) this imphes W^y = W^y. A further 
restriction is provided by the current conservation q^W^y = q^W^y = 0. 
The general expression satisfying these conditions can be written as 

where Wi and W2 are, a priori unknown coefficients. It is traditional to 
define dimensionless structure functions 

Fi(x,g2) = MTVi F2{x,Q^) = vW2, (2.7) 
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which, in the <^ hmit, can be projected from the hadronic tensor as 



F2 

X 

Fi 



1 2t2 

-g^"" + -^^P^P"" 1 MW, 



(2.8) 



1 



In terms of the structure functions Fi and F2 the cross-section in Eq. fl2.3p 
can expressed in an invariant way 



dxdQ"^ 



47rc4^1 

rc 



xy^Fi + F2 1 - y 



xyM'^ 



where s = {P + k)"^ denotes the center-of-mass energy, and 
stands for the fine-structure constant. 



(2.9) 



eV47r 



Parton model 

The parton model [H [5] can be motivated by considering the DIS not in the 
target-rest-frame but in the electron-proton center-of-mass system. In such 
a frame, the nucleon appears Lorentz contracted, and the time dilatation 
slows down the intrinsic interaction rate of the fundamental constituents of 
the nucleon, the partons. During the short period it takes for the electron to 
traverse across the nucleon, the state of the nucleon wave function can thus 
be envisioned as being frozen to a superposition of free partons colhnear with 
the nucleon. Mathematically, the parton model is defined by the relation 

da = y] C d^dalimm). (2.10) 

q Jo 

where o-q(^P) is the leading order (Born) cross-section for the electron- 
parton scattering, with the parton carrying a momentum p = ^P. The 
functions fq{^) are called parton distributions, and represent the number 
density of partons of flavor q in the nucleon. 

7 q 



Figure 2.2: The leading-order diagram for photon-quark interaction. 
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In QCD, only quarks carry an electric charge Cq and the definition ()2.10p 
with Eq. fl2.3p imphes that the hadronic tensor W^jy can be written as 

^f^'^ = E r T^wo, (2.11) 

where its partonic counterpart W^u is essentially the square of the diagram 
in Fig. ^ 

= f-Q^Tr[/777l<^(e-x). (2.12) 

Neglecting the nucleon mass M compared to the photon virtuahty Q^, 

g^, TvyrpY] = P^P. Tr[/7^^7n = 0, (2.13) 

we find 

- g^^AMW^^,) = elx5{^ - x) pf^p^MW^,) = 0. (2.14) 

Consequently, the parton model predictions for the structure functions re- 
duce to an electric-charge-weighted sum of the quark distributions, 

2xFi{x) = F2{x) = Y,^lxfq{x), (2.15) 
q 

and the cross-section in Eq. (12. 9 p can be written as 
where (Tq denotes the partonic Born cross-section 



2 V " S-M2 



(2.17) 



It is a prediction of the parton model that the structure functions Fi 2 are 
only functions of x, and should not depend on in the » hmit. 
This phenomenon, termed as Bjorken- scaling, was indeed observed in the 
early SLAC experiments providing direct evidence about the inner con- 
stituents of the nucleon. Later experiments which have covered a larger do- 
main in the {x, Q^)-p\ane have revealed, however, that the Q^-independence 
of the structure functions Fi 2, although a good first approximation, was not 



8 



exactly true. Such deviations are clear e.g. in Fig. 12. 3[ which shows some 
experimental data for the proton structure function F2. The scaling viola- 
tions, as they are nowadays called, can however be fully explained by the 
QCD dynamics — by so-called DGLAP equations, to be derived shortly. 
Together with the asymptotic freedom, these equations with factorization 
theorem constitute the main pillars of perturbative QCD. 




Figure 2.3: Experimental data for proton structure function F2 from SLAG 
^ and CERN-EMC experiments. 



2.2 Initial state radiation 
2.2.1 Origin of the scaling violations 

Due to the inclusive nature of the deeply inelastic scattering nothing forbids 
having additional QCD particles in the final state. First such corrections to 
the Born-level matrix element originate from a radiation of a real gluon as 
shown in Fig. 12.41 Both of these diagrams are divergent as the intermediate 



g 




p. J 

Figure 2.4: Real gluon radiation. The letter i,j,a are the color indices. 

quark propagators are close to being on-shell: 

{p-k f = -2/fc°(l -cos^) ^0, 
[p' + kf = -2p'°k° (1 - cos 0') 0. 

This can happen either if the momentum of the emitted particle goes to 
zero, k^,p'^ 0, or if the emission is in the direction of the incoming or the 
outgoing quark 9, 6' 0. These are archetypes of infrared and collinear 
singularities correspondingly. There are also same kind of divergences stem- 
ming from the virtual corrections, and it turns out that all but the collinear 
divergence related to the gluon radiation from the incoming quark will 
eventually cancel. In what follows, I will show how to extract these diver- 
gences and how their resummation gives rise to the parton DGLAP evolution 
— the Q^-dependence of the parton distributions observed in the experi- 
ments. 

2.2.2 One gluon emission 

Rather than drawing graphs for matrix elements, for the rest of this Chapter, 
1 will draw the graphs in the cut diagram notation (see e.g. [8j) directly for 
the cross-sections. The square of the diagrams in Fig. 12.41 look as 




Figure 2.5: Diagrams representing the 7*5 qg process. It should be 
understood that the parton lines in the middle are real, on-shell particles. 

Although the such squared matrix element is certainly gauge invariant, the 
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contribution of an individual graph depends on the choice of gauge. How- 
ever, as I already mentioned it is only the collinear (p — k)'^ ^ singularity 
which turns out to be the relevant one. In the Feynman gauge, all but the 
second of the diagrams in Fig. 12.51 will contribute in this kinematical hmit, 
but, as it turns out, in the axial gauge it is the first diagram alone that 
is responsible for the divergent behaviour. For obvious reasons I call it a 
ladder diagram. 



Axial gauge 

The class of axial gauges is specified by a gauge- fixing term — ^ (n • G)^ 
in the QCD Lagrangian where G denotes the gluon field, n is an arbitrary 
four-vector and ^ is the gauge parameter. The gluon propagator in this 
gauge is 



D^u{k) — ^2 



K ■ n [k ■ nj^ 



(2.18) 



The sum over the two physical polarization states e^^''^{k) (k"^ = 0), obeying 

k ■ €^^{k) = and n ■ e^^{k) = 0, normalized by e^^{k) ■ e^^{k) = —1, reads 

A 

Usually, it is convenient to choose ^ = and re^ = which specifies the light- 
cone gauge. The axial gauges are sometimes called physical gauges: the 
reason for this is most distinct in the hght-cone gauge as any propagator 
in a Feynman diagram can be replaced by the polarization sum over the 
physical states: 

A convenient choice for the light-like axial vector in the present problem is 

n = q + r]p, r] = . (2.21) 

2p-q 



Sudakov decomposition 

In extracting the dominant part of the squared matrix elements, it is con- 
venient to parametrize the momenta of the outgoing partons by [9| 

k = {l-z)p + Pn + k±, (3 = — ^ , (2.22) 

2(1 — z)p ■ q 
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where k± is a space-like 4-vector orthogonal to n and p: < 0, n ■ kx = 
p ■ kx = 0. For example, in the center-of-mass frame of p and n, 

p = (P,0,0,P) 

n = (P,0,0,-P) (2.23) 
k± = (0,k^,0) 

where P is some reference momentum. In such a frame the interpretation 
of /e^ = — k^^ as the transverse momentum is evident. Furthermore, 

t' = ip-k)' = kl/{l-z), 

and we see that the colhnear t ^ divergence should be found by extracting 
the l/Zc^-pole. 

The calculation 

The squared matrix element corresponding to the first diagram in Fig. 12.51 
reads 

l>^^;*1i,^J^. = C^F5s'|EF^[^V/K/7.], t = p-k (2.24) 

pol 

where the color factor Ci? = 4/3 arises from (see Fig. 12. 4p 

i,j,a a 

Using the polarization sum Eq. p.l9p one finds 

^^l^^ = -l-(^ + /3^), (2.25) 

pol 

and after a short calculation 



/(^ + /3yi)/ = 




(2.26) 



where the remaining terms are higher order in k± and will not contribute 
to the colhnear divergence. In total, 

\M};t^%\l. = dCF^-^^ (\^) X |tV [/7,|^7.] + • • • • (2.27) 

12 



It is essential that the last combination of terms is nothing but the squared 
matrix element in the Born approximation. Supplying the phase-space ele- 
ment in the Sudakov variables 

(fk 1 dz „ 2 



(27r)32A;0 lOvr^ 1 - z 

one obtains 



dk±, (2.28) 



'^'^ n-^7*™l^. = #^-(S)^^^(^)><|Tr[/7,^7.]+---, (2.29) 



(27r)32A;0' ^2 V27r. 

where 

(2.30) 



1 + z^ 



is the so-called Altarelli-Parisi splitting function associated with the un- 
polarized quark quark transition. In the colhnear hmit, the variable z 
is readily interpreted as the momentum fraction of the quark left after the 
gluon emission. The contribution to the quark tensor Wjiu is 



Neglecting all 0{k'j_) terms which would cancel the collinear singularity in 
Eq. ([CT]) . 

p" = 2z{p.q)+q' + 0{kl)^Q' 



X 



and 

g^, Try'j^pY] ~ -4Q'| P^Pu Try'j^p-f"] = 0. (2.31) 
Thus, the dominant 0{as) piece in the quark tensor is 



+ (2.32) 



which contributes to the hadronic tensor by 

'Q,\ f dk^ dz 

9 



(2.33) 



As anticipated, the collinear divergence manifests itself in the / dk\ /k^ in- 
tegral. The upper bound for this integral is proportional to but the lower 
hmit remains zero for massless quarks. Even if the quark had a small reg- 
ulating mass m, the resulting logarithm log((5^/m^) would not be infrared 
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safe: the resulting cross-section would be sensitive to the value of in 
the large-Q^ limit. The solution to this problem will require resumming a 
whole tower of such logarithms. For the time being, however, I add this 
divergent piece to the DIS cross-section 



d^a LL d (To 



dxdQ^ 



dxdQ"^ ^ 



( OLs 

1 + — 
V27r 



1 I 

log — 



(2.34) 



where the designation LL means that 1 have kept only the leading logarith- 
mic contribution, and the shorthand notation ® stands for the convolution 



^ dz 
—6{1 

z 



qq: 



fq{x) 



(2.35) 



Since the left-hand side of Eq. (I2.34p is a measurable, finite, quantity the 
non-perturbative parton density fq is inevitably intertwined with the arbi- 
trary cut-off scale such that the cross-section is finite. 



1 still need to prove my claim that in the axial gauge this is the only colhnear 
logarithm related to the initial state gluon radiation. This is actually quite 
a simple task: writing down the cross term 




pol 



one reahzes that if the trace is independent of k± there will be a similar 
colhnear logarithm as found above. However, noting that 



--2z{p-e), 

the polarization sum reveals the structure 

demonstrating that no /c^-independent term exists and the proof is com- 
plet£l. Thus, in the colhnear hmit and in the axial gauge, there is no 



^In the Feynman gauge with e^(fe)ej^(fc) = —g^v this last step would not be true. 
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interference with the outgoing quark and, in the spirit of parton model, the 
factor 



can be interpreted as a probabihty density for the quark to radiate a gluon 
carrying a fraction 1 — 2; of the quark momentum, before getting struck by 
the photon. One should note that due to the 1/(1 — 2;)-pole in Pqq{z), this 
probability diverges in the infrared z ^ 1 hmit, making the convolution 
integrals apparently ill-defined. However, the probability that the quark 
re-absorbs the emitted gluon diverges similarly and will wash out the z = 1 
singularity, as will be discussed soon. 

2.2.3 Multiple gluon emissions 

Based on the previous section, it is natural to expect to find two similar 
colhnear divergences as in Eq. (12.341) if double gluon emission, shown in 
Fig. 12. 6[ is considered. This is indeed the case and employing the method 
introduced earher one can extract an log^((5^/m^) contribution to cross- 
section Eq. (I2.34p . This is how it goes. 



Figure 2.6: Ladder graph for two-gluon emission. Momenta of the produced 
gluons with polarizations £1,62 are denoted by ki,k2, and the intermediate 
quark momenta are ti = p — ki, t2 = p — ki — k2. The color indices are 
denoted by i,j,h,ai,a2. 

The squared and spin-summed matrix element for the ladder diagram in 
Fig. [22] reads 



where the color factor Cp = (4/3)^ arises in the following way (summation 







pol 
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over all indices is implicit): 

1 



lxr(t«it'^i)iTr(t"2t«2) = 



(2.38) 



r2 



where I used {t"-t°-)hs = (1/3) ^/i^. Introducing the Sudakov decomposition 
for the lower gluon momentum 



(1 - zi)p + Pin + k 



Pi 



2(1 - zi)p -q' 

one immediately obtains, reading from the preceding calculation, that 

poll 



(2.39) 



1 + zf 
Zl \ l- Zi 



{-kli_)i> + . . 



(2.40) 



where 1 have again omitted the terms higher order in fci^. In the same way, 
writing the Sudakov decomposition for the upper gluon momentum as 



k2 = Zi{l- Z2)p + I32n + fc2±, 



/32 



2zi{l - Z2)p ■ q 



(2.41) 



and dropping terms higher order in ki± and /c2±, one finds the leading 
contribution 



P0I2 



l + zl 



1 - Z2 V 1 - -^2 



(2.42) 



Thus, the squared matrix element (12.370 acquires a form 



\M 



Ladder 1 2 
'y*q^q,2g\fiu 



ys ,4 



1-zi 



2P,,(^2) 
1-Z2 



(2.43) 



where the last factor is again the Born matrix-element that has penetrated 
through the calculation. If there were not the factors if 2 the denominator, 
the leading factors for both emitted gluons would be identical. However, 



l-zi 

^2± 
1-Z2 



1-^1(1-2:2), 2 , I, 1 

H 7-2 k(± + h± ■ k2±, 

1 Z\ 
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(2.44) 



where the latter one looks bad. In the region of phase space where —k 



11. 



< 



-/cg^ one can power expand Eq. fl2.43p in k\^/k2 



2±i 



schematically 



1 1 



I A^Ladder |2 

\J'^l-^*q^q,2g\iJLy ^ ,2 7,2 



'^2± 



+ B 



A-2 



+ 



(2.45) 



where the odd powers of /cij, are absent as they would vanish upon integra- 
tion. Whereas the integration over the first term gives the leading double 
logarithm, 



i2 



1 

2! 



log' 



(2.46) 



the rest can give only a single logarithm. In the opposite transverse mo- 
mentum ordering —klj_ > —k^^, one again obtains only single logarithms. 
Thus, the leading contribution stems from the transverse momentum order- 




Figure 2.7: Diagrams for two-gluon emission that do not contain double- 
logarithms. 

ing —kfj_ ^ — emitted gluons, and 



I A^Ladder |2 LL — ^1 ^2 -n / . \1 ~ ^2 



[29lPqq{z,)\ {2glPqq[z2)\ (2.47) 



Following the same steps as earlier, we find 



1 



« 2 



zvr V 



2=122 



dz2 

Z2 



dzi 
z\ 



(2.48) 



which contributes to the hadronic tensor by 



j:4 



zvr 
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-Pgg ® A<? ® /g- (2-49) 



The convolution between three objects above is defined by 

Pqq®Pqq®U=f —Pgq{z2) [ —Pgq{zi)f, 
Jx ^2 Jx/zo ^1 



(2.50) 



with obvious extension to convolutions between an arbitrary number of 
functions. Thus, to 0{a1), the leading logarithms organize themselves as 



d cr LL d 170 



dxdQ^ dxdQ' 



2 \2irl * 



^]P 



(2.51) 



Pr 



qq 



fq- 



Based on a similar reasoning as in the end of the previous subsection, the 
diagrams like those in Fig. 12.71 cannot contain ©(a^ log^(Q^/m^)) terms in 
the axial gauge — it is the ladder diagram in Fig. 12.61 alone that gives the 
leading logarithmic singularity. 

The generahzation to an arbitrary number of colhnear gluon emissions from 
the initial quark is now quite straightforward: For n emitted gluons the 
leading logarithms originate from the region of the phase space where the 
transverse momenta are strongly ordered 

and the contribution to the DIS cross-section is 




k„ 1 



Figure 2.8: Ladder graph for n-gluon emission. 



^ E m" -og" ( « ■ ■ ■ P,, ^U (2.52) 



n times 



18 



Thus, the leading logarithm contributions to the DIS cross-section constitute 
a series which is formally an exponential 



dxdQ^ dxdQ' 



E 



exp 



(2.53) 



Comparing this expression to the corresponding parton model prediction, 
given in Eq. (j2.16p , one can see that the resummation of the leading log- 
arithms is equivalent to replacing the Q^-independent parton distribution 
function by 



fq{x) fq{x,Q^) = exp 



27r 



log 



9 I ^ qq 



fq 



(2.54) 



Taking the Q^-derivative we see that Q"^) satisfies the following integro- 
differential equation 



d 



Pqq'^fq{x,Q^ 



(2.55) 



which is an archetype of the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
evolution equations [lOl dH SSI [13] , or DGLAP equations in brief. 



2.2.4 More splitting functions 

The gluon emission discussed above is, of course, only one possibility among 
other QCD-interactions. For example, from 0{as) onwards, also gluon- 
initiated subprocesses contribute to the deeply inelastic cross-section. The 
simplest such diagram is shown in Fig. 12.91 Similar to the gluon radiation 




Figure 2.9: A gluon-initiated ladder diagram. 

graphs discussed in the preceding sections, also this diagram — and in the 
axial gauge this ladder-type diagram only — gives a collinear divergence. 
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Extracting this divergence is rather straightforward having already care- 
fully done the groundworks during the previous sections. Starting from the 
squared matrix element 

\M};t^%-Xu=9'sTRellY.Ji^[j^'^^^^^^^^^^ t = P-k (2.56) 

pol 

where 

is the appropriate color factor, and introducing the Sudakov decomposition 
(I2.22I) for the outgoing quark momentum k, the polarization sum gives 

^ = 2 [(1 - z)^ + . (2.57) 

pol 

To leading power in A;^ we find 

/ [(1 - z)p + Pi]i = ^ [(1 - zf + z^] p, (2.58) 
and consequently 

\M};t^%-Xu = 9s^^^ Tn [(1 - zf + z'] X |Tr [/7,^7.] • (2-59) 

Comparing to Eq. (12.271) we reahze that the leading contribution to the 
deeply inelastic cross-section from the graph 12.91 becomes 




(2.60) 



where 

P,g{z)=TR[z^ + {l-zf] (2.61) 

is the splitting function for a gluon^quark transition and fg is the parton 
distribution function for the gluons. 

Having now considered two different ladder vertices, we can also pile them 
on top of each other to form a parton ladder like the one in Fig. l2.10l below. 
The corresponding squared matrix element reads 

\M^*G%-gg\lu = 9sTrcA J2 ilV [^'^MtWi^MtM ' (2-62) 

^ pol ^1^2 
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f2 = p - fcl - k2,j 

ti = p — ki, i 



T5T5T5T5T5T5T5T5Tr 



Figure 2.10: Another gluon-initiated ladder diagram. 



where the color factor arises as 

^ ' 4.b j.a \ /4.b j.a \* 



= -Tr(t'^t'*)-Tr(t¥) 
8 3 



Using the momentum decomposition ()2.39p , the inner part of the trace above 
gives the leading factor 



E 

poll 



2— [{l-zif + zi]f+..., 



(2.63) 



and what is left inside the trace is, using ()2.4ip , 



pol2 



i-kl^) (l + zl 



I - Z2 \l - Z2 



+ .... 



(2.64) 



The leading region is again the one with —klj_ < —k"^^, and therefore 



\M 



Ladder |2 LL 1 — 2:2 



2 [2glPqq{z2)\ WsPqgizi)] (2.65) 



-k 

|Tr [/ 



i± 



which leads to a contribution 

d'^f^Q sr^ o 1 / a., \ 2 



dxdQ'^ 
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to the deeply inelastic cross-section. 



Thus far we have only considered parton ladders with quarks at vertical 
hnes. To illustrate how the calculation works when there are vertical gluon 
Unes, let us consider the diagram shown in Fig. 12.111 The squared matrix 




t2=p-ki - k2,j 



ti = p — ki.a 




Figure 2.11: Ladder diagram with gluon as a vertical hne. 



element is in this case 



\M 



Ladder 1 2 



e 1 

glCpTR^ ^ 7474 \^lii^hoM2l^tklv\ Tr [^i 

^ pol ^1^2 



t\ • n 



t\ ■ n 



(2.66) 



where the color factor comes from 



1 



3 

-Tr(t"t")-Tr(t¥) 
3 8 



Applying the Sudakov decomposition (12.391) , we find 



8 ^ -kl^ 
zi '-^ 1 - zi 



ti ■ n 



5"" + 



ti ■ n 



p ■ n 



+ • 



(2.67) 



Whereas the second term in the lower line explicitly contains the sum over 
the vertical gluon polarization states, the first term looks a bit puzzhng. The 
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trick is to notice that under the integration over the transverse momentum 
this term becomes 



/ 



rtjj^ _|_ 



p ■ n 



D{kl 



(2.68) 



where D is a function that depends only on k\ and not on the transverse 
components separately. Thus, we may write 



ti ■ n 



ti ■ n 



p ■ n 



(2.69) 



Applying again the Sudakov decomposition ()2.4ip we have 



tklaH2lr]tk -5"" + 



p ■ n 



"2^ 



1-Z2 



[(1-^2)^ + ^2] zip+.^ 



and in the leading region 



Ladder 1 2 
'y*q' ~*q'qq I fiu 



- [29!Pq,{z2)] ^ [29!P,q{z,)] (2.70) 



-k 



'21. 



where 



l + {l-zf 



Pgq{z) = Cp 



(2.71) 



is the splitting function for the quark^gluon transition. The squared matrix 
element above now leads to a term 



in the deeply inelastic cross-section. 



The remaining splitting function to be calculated is Pgg corresponding to 
the gluon^gluon transition. This can be computed from the ladder diagram 
depicted in Fig. 12.121 which corresponds to the squared matrix element 

„2 



pol 



(2.72) 



l''2 



/"(ti + pf + + k^f + /^(ti - k^y^ 

t{nP + 1{ 



ti • n 



ti ■ n 
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t2 = P — kl — k'J: j 



t\ = p — kl, a 



Figure 2.12: A ladder diagram from which one can compute sthe plitting 
function for gluon^gluon transition. 



The evaluation of this matrix element is somewhat more tedious than the 
previous ones, yet straighforward. The resulting leading logarithmic contri- 
bution to the deeply inelastic cross-section is 



with 



Pgg{z) = 2Cg 



1 — z z 

+ - + z{l-z) 

z I — z 



(2.73) 



Prom now on, one can pretty much see how this goes on: each additional 
ladder-compartment in which parton of flavor i transforms to j, effectively 
just increments the power of OslogQ^/m^ by one unit and adds the cor- 
responding splitting function Pij to the convolution integral. The possi- 
ble building blocks for constructing the ladders are displayed in Fig. 12.131 
together with the characteristic splitting functions. Clearly, we should ac- 
count for all possible parton ladders — also the gluon-triggered ones — 
when defining the scale-dependent quark densities. Therefore, we define 
the scale-dependent parton distributions as a sum of all possible ladders 
that end up with the specific parton. The appropriate generahzation of the 
definition ()2.54p can be neatly written down as a matrix equation 

( ) , (2.74) 

where qj should be understood as being a vector with different quark fiavors 
as its components and the sphtting functions Pq^qj, Pq^g, Pgqj as matrices 



fa{x,Q'') 



exp 



as, Q' 



Pr 



P, 



P. 
P 



.9.9 
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Cp = 4/3, Ta = 1/2= = 3 




Figure 2.13: Unpolarized splitting functions. 

with the appropriate dimension. In the leading logarithm approximation, 
the sphtting functions are fiavor-bhnd and we can write, exphcitly 



P P 



99] 



99 



( p,. 





• 


■ P,9 \ 





p 

-Tqq 


• 


P 








p 

J^qq 


p 

^qg 



V Pgq 



p p 
^9<i ^gq 



^99 / 



(2.75) 



The complete set of DGLAP evolution equations follow by taking the Q"^- 
derivative: 



as 
27r 



P P 

^qiqj ^qi9 



gqj 
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fqAQ') 

fgiQ') 



(2.76) 



In summary, the leading collinear singularities in the perturbative Feynman- 
diagram expansion can be factored to the scale-dependent parton distribu- 
tions fiixjQ"^) such that the parton model prediction for the DIS cross- 
section stays formally intact, but the parton densities no longer respect the 
Bjorken-scahng but are Q^-dependent. Also the interpretation of the parton 
distributions as simple number densities upgrades to being number densi- 
ties with transverse momentum up to Q^. This and further extensions to 
the simple parton model are often referred to as pQCD-improved parton 
model. 

There are still, however, two serious deficiencies in the equations ()2.76p : 

• The sphtting functions Pqq{z) and Pgg{z) diverge as z ^ 1 due to the 
l/(l-2;)-poles, making the convolution integrals thereby meaningless. 

• The argument of the strong coupling constant as remains undefined. 

In order to fill these gaps, we need to discuss also the virtual corrections to 
the parton ladder on a same footing with the real parton radiation. 



2.3 Virtual corrections 
2.3.1 Quark self-energy 

The Sudakov decomposition of the momentum provided a useful tool for 
extracting the colhnear hmits of the parton radiation diagrams. This is 
also true in the case of loop-integrations, but the parametrization must be 
shghtly modified to account for the non-zero virtuality of the loop momen- 
tum k. Also, the virtuahty of the parton (as it may lie in the middle of 
the ladder) should be kept arbitrary. As an appropriate extension to ()2.22p , 
1 will decompose the loop momentum a^ 

k = zp + Pn + k±, p = — ^ — — ^. (2.77) 

2z[p • n) 

The axial-gauge expression for the quark self-energy diagram shown in 
Fig. [2lH reads 



(27r)4 [[p - ky + ieWk"^ + ie] 



■ „^ k^'n'' + k''n^' 

9 Z 

k ■ n 



^Fcynman(p) ~l~ ^Axial {?)■ (2.78) 



^An explicit definition of the axial vector n is not needed here. 
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p p-k p 



Figure 2.14: Quark self-energy correction at 1-loop. The momentum of the 
parent quark is p while the loop momentum is denoted by k. The letters 
i,j,a stand for the colors indices. 



I will now carefully show how to extract the colhnear divergence starting 
with the Feynman gauge contribution 



^Feynman(p) — 

First, in the numerator 



(27r)4 [(p- A:)2 +ie][A;2 +ie] 



(2.79) 



-2{f - 

-2 [{l-z)f- 
-2 [{l-z)f- 



(2.80) 



where I dropped the ^j^-term as its contribution will vanish as an odd inte- 
gral. The denominator from the quark propagator is 



{p — k)'^ + ie 
where 1 defined 



1 - z 



k' 



k^, 



1 - z 



zp^ — ie' 



1 - z 



<0if2;<0or2;>l 
> if < z < 1. 



(2.81) 



(2.82) 



The sign of e' is essential in defining the locations of poles in the A;^-plane. 

z 1 



SFeynman(p) = - ^Cp^ / dzdk\dk'^^ — 
ZTT ZTT J \Z\ ^ 

(1 - z)p - I3i 



[A:2 + 



(2.83) 



Except for the /c^-terms in the numerator which will not contribute to the 
A;^ ^ divergence, the fe^-integral can be evaluated as a contour integral. 
Depending on the value of z, the locations of the fe^-poles and the convenient 
integration contours are shown in Fig. ()2.15p . If < z < 1, the integration 
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Figure 2.15: The integration contours and the pole structure for < z < 1 
(left) and z <0, z > 1 (right). 



contour can be closed in the lower half-plane enclosing k'^ 



-ie pole 



oo [k"^ + ie] k"^ — — zp"^ — ie' 



-27ri 



l-z 



k2 



_i_ — z{l — z)p'^ ' 



(2.84) 



while for other values of z the integral contour can be closed in the upper 
half-plane where there are no poles and integral gives zero. Therefore, 



eynman 



ip) 



V27r 



I \ — ) Cp dz 



k5_ - z{l - z)p'^ 



(2.85) 



as 
2tt 



Cf log 



2z{p ■ n 

A2 ■ 



dz{l — z), 



where 1 have dropped all but the logarithmic contribution diverging as 

0, and regulated the UV-divergence by a cut-off A^. The axial contribution 



d'^k 



+ 



(27r)4 [{p - kf + ie][A;2 + ie] k ■ n 



(2.86) 



can dealt with similarly. Applying the Sudakov decomposition plus drop- 
ping all fc^-terms and those ones odd in k±, the numerator simphfies to 



+ 



(2.87) 



Apart from the factor k ■ n = zp ■ n, the denominator is identical with the 
Feynman-gauge case and the fc^-integral can be performed in a similar way. 
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The result is 

S Axial (p) 



p ■ n V 27r 



k2 



2p ■ n V 27r 



CFlog 



- z{l — z)p'^ 
dz- 



(2.88) 







1 



z 



Thus, the total self-energy correction is of the form 



T,{p) = Sf, 



cynman 



(P) + ^Axialb) 



2p ■ n 



ifi 



and the 1-loop corrected quark propagator becomes 

pi pi pA pz zp ■ n 



i. 



(2.89) 



The last term above will lead to a loss of a leading logarithm in the par- 
ton ladder and thus the 1-loop correction to the quark propagator can be 
effectively accounted by a multiphcative renormahzation constant 



\- A-B 



OLs_ 

27r 



A2_ 



dz 



1 + z' 



(2.90) 



This is also the residue of the ^ propagator pole needed (through 
the LSZ-reduction) if the quark is an on-shell final/initial state particle. 
Although 1 did not carefully keep track of the UV-divergences, the Eq. (12.901) 
is also accurate in this sense — it is a cut-off regulated version of the result 
regulated by going to 4 — e dimension^ 



1 



V27r 



+ log 



Cf / dz 



1 + z' 
1-z 



e v-P 

To understand how this damps the singular behaviour of Pqq{z), it is sim- 
plest to look at the effect of applying the 0{as) external leg correction to 
the leading order contribution in the DIS cross-section (I2.34p . 

.27r. 



da 



LL 



ddo ^ ( 
ddo ^ ( 



Pr 



In, ( 



P 



<?<? 



where I defined the regulated sphtting function Pqq{z) by 



Pqq{z) 



Cf 



1 + z' 



l + z 



1 - z 

2 



z 



6{1 



dx 







(2.91) 

(2.92) 
(2.93) 



^The loop calculations in the light-cone gauge are also little tricky, see [14[ [TF 
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Here, we meet so-called plus distribution, which should be understood 
through integration against sufficiently smooth "test-function" h{x): 



dzh{z) 



dz 



l + z' 



[h{z) - h{l)] 



(2.94) 



There are two important things to be emphasized: First, the inclusion of 
the virtual correction serves to regulate the z ^ 1 singularity of the spht- 
ting function Pqq{z), as promised. Second, the virtual piece m^) gets 
replaced by Zq{—Q'^) which no longer contains colhnear ^ divergence. 



2.3.2 Gluon self-energy 




Figure 2.16: Diagrams contributing to the gluon self-energy correction at 
1-loop in the axial gauge. 

The gluon-self energy 1-loop correction can be calculated following essen- 
tially the same procedure. In the axial gauge, there are only two contributing 
diagrams, shown in Fig. ()2.16p . Evaluating the quark and gluon loops the 
following logarithmic pieces are found 



•«s 2 , 



P^iPu \ 1 / AM 2 



log ( ] Cg 



27r 



2 / P^lPv \ Al „ 



+ 



-n 



p ■ n 

Thus, the result is of the form 



p. 



W _ J^n"" 12(1 



P 



p ■ n 



(2.95) 
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AP 9iiiy 



P^iPu 
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p ■ n 



P' u 

n 



p ■ n 



inducing a correction 



A;^ \ k ■ n 



ki'ni + k1n^' 
k ■ n 



iB 



{k ■ n) 



k^n" + k^n^ 
k ■ n 
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to the gluon propagator. Again, the latter part above does not give a leading 
logarithm, and the loop insertions can again be effectively accounted for by 
a multiplicative renormalization constant 



1 , 1 



A2_ 

_p5 



dz 



2 T ■ 



(2.96) 



As a consequence, the gluon gluon sphtting function Pgg{z) gets replaced 
by a regulated one 



Pggiz)=2CG 



1 



z {l-z)+ 



+ z(l 



+ 



11 

"6 



■Cg - i^rifTj 



6{l-z). 
(2.97) 



2.3.3 Renormalization of the ladder vertex 

The standard field theory text books (e.g. [TBlCLTj) relate the running cou- 
pling constant (7s (^^) to the bare one gos e.g. by 



with k\ = k"^ = = — and where Zi is the renormalization factor for 
the qqg-vertex. At 1-loop, we have the well-known result 



Q,0 

47r 



a2 . , 



/3olog /xVA? 



(2.99) 



where /3o = ^Cg — ^T^nf, A^ in the first hne denotes the cut-off, and 
Aqcd ~ 200 MeV is the QCD scale parameter. 



In a parton ladder, however, the kinematics appear quite different from that 
in (12.981) — the virtualities k\,k\ of the partons in the vertical hues are 

,2 
C 



strongly ordered while the horizontal rung k"^ is on-shell: 



- k% -k% > 0, -k% < -k%, kfj = 0- (2.100) 

In order to discuss what is involved, 1 will adopt a specific example with 
quark (A) sphtting to a on-shell quark (C) and virtual gluon (B). The out- 
going quark contributes by a factor Zq{—m?) with m? 0. However, to 
the same order in Oia.) we should also consider the contribution from the 
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Figure 2.17: Various renormalization factors for a ladder vertex. 

o 



Figure 2.18: Diagram for gluon radiation from the outgoing quark. 

horizontal quark rung radiating an additional gluon as shown in Fig. 12.181 
A calculation along the hues presented in Sec. I2.2.2[ reveals that the leading 
contribution is effectively a multiplicative factor 

[Ladder skeleton] x Cp— / -r^ / dz , (2.101) 

2vr Jm^ k|_L Jo 1- z 

where the upper limit k^ denotes the aggregate transverse momentum car- 
ried by the outgoing quark-gluon system and derives from the requirement 
not to upset the underlying logarithmic structure of the ladder. This is 
exactly of the same form as the quark renormahzation factor Zq{—m?), and 
when combined, the sum is clearly free from colhnear divergences. In ef- 
fect, we may simply make a replacement Zq{—m'^) Zq{k'^). Similarly, 
the incoming quark line involves a renormahzation factor Zq{k\) which gets 
replaced by Zq{k'^) by the mechanism demonstrated in Eq. ()2.9ip when the 
contribution of gluon radiation is included. 

The remaining piece is the vertex part Zi{k\,k'^,k(^). However, in the 
axial gauge it happens that if k"^ is kept finite, Zi does not contain terms 
that would be divergent in the /c^ ^ ^ limit. In other words, in axial 
gauge all mass singularities are contained in the self-energy factors and we 
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may safely replace Zi{k\,k'^,k'^) by Zi{k'^,k'^,k'^) without losing large 
logarithms. Thus, the virtual corrections indeed eventually combine to the 
usual definition of the running couphng, 



Zq{ky] 



z,{kl) 



ZG{kl)Jz,{kl] 



Zq{ky] 



Zi{k\, fc^, k'i.) 



Zi{k'^, k'^, k'^j 



9s{-kl). 



Incorporation of the running coupling to the resummation of leading loga- 
rithms is straightforward: In each ladder vertex we change as(k^), 
and do the nested transverse momentum integrals like (12.461) by a change of 
variables 



Pi 



log 



a.(ki) 



(2.102) 



such that 



,2 k^^ 27r 



f: 



2x dk{^ as(kj^) 



k2 



27r 



(2.103) 



The correspondingly corrected definition for the scale-dependent parton 
densities becomes 



exp 



Pr 



HQj 



9<lj 



P 



99 



(2.104) 



where the sphtting functions are understood as being the regulated ones. 
By differentiating with respect to k and changing variables, we obtain the 
complete DGLAP evolution equations that resum all leading logarithms: 



Q 



Q 



5Q2 



2dfg{x,Q^) 



27r 



[Pg.g,^fgAQ') + P9.9^f9iQ' 



(2.105) 



2.4 Higher orders 



In this section I have systematically kept to the leading logarithm approxi- 
mation, in which the logic is to retain only terms of the form a" log"((5^/m^), 
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discarding all contributions which are suppressed by additional powers of 
Ug. However, if one keeps track also of the non- leading contributions 

log"(QV"^'), log"(QV"i'), • • • 

one finds that the splitting functions Pf^f^ actually constitute a power series 
in as, 

Of these, P(^) have been known since 1980's [IH [l9] but P^^) jj^ve been 
computed only recently [20l[2l]. While the kernels pf^ are unique, inde- 
pendent of how the colhnear divergences are regulated, the higher order 

in) 

splitting functions Plj , n > are not unique but they depend on the 
framework. 



2.5 Factorization in deeply inelastic scattering 

For couple of times during this section, I noted that when the divergent 
logarithms from the parton ladders were extracted, a multiplicative term 
that coincided with the leading order cross-section for photon-quark scat- 
tering, was found. This property is not, however, a special feature of the 
leading order cross-section, but order by order in perturbative calculations, 
the divergent logarithms can be systematically factored apart from the per- 
turbative parton-level pieces 



da. = daT + (^j dat^ + ^^j da^' + • . . (2.107) 

That is, the cross-section for deeply inelastic scattering retains its same 
simple form 

i 

where the parton densities fqiQ"^) obey the DGLAP equations ()2.106p . This 
is a special case of the pQCD factorization theorem [22l [23] which states 
that the colhnear singularities can be, order by order, factored to the scale 
dependent parton densities, and finiteness of (I2.108P is guaranteed. A more 
complete treatment indicates that the factorization is subject to power cor- 
rections O (^^qcd/Q^^ which may become important for small Q^. Such 
terms appear, for example, if the partons are not exactly collinear with 
the parent nucleon, but are allowed to carry some "primordial" transverse 
momentum k^. More generally, such terms arise form multi-parton inter- 
actions. 
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The definition of the parton densities is not unique: starting from parton 
densities fi{x,Q^) we may define another version f[{x,Q'^) by 

/;(x,Q2)^^C,,^/,(Q2), (2.109) 



where C,j{z) = 1 + ^c\]\z) + {^f c\f{z) + . . .. In terms of the primed 
densities, the cross-section can be written as 

where I defined cj,' = Cij^ ■ The same reshuffling imphes that obey 
the DGLAP equations with sphtting functions 



An 

C®P®C-^ - 27r/?(a,)— 



In other words, beyond leading order, the perturbative coefficients o-j, and 
the splitting functions P are intertwined with the exact definition for the 
parton densities. These kind of different "bases" for computing are known as 
factorization schemes and, in principle, the predictions for physical, mea- 
surable cross-sections should be independent of the chosen scheme. How- 
ever, as the change from a scheme to another is facilitated by the perturba- 
tive coefficients Cij — often motivated by some particular cross-sections — 
the predictions from various schemes are not precisely equal. However, such 
difference is always one power higher in Ug than to which the computation 
was performed. 



There is a similar ambiguity in choosing the argument of /g(Q^) in ()2.108p . 
This is because different scales are related perturbatively as 

h{x^Q^) = h{x,Q)) + ^\og{Q''/Q))Pi,®f,{Q)) (2.110) 



+ 



2 V27r/ 



Y.Dij{Q''/Q))®f^{Q)), (2.111) 



and defining &j{x,Q'^ jQj-) = d-i{Q'^,Q'j) (S!) Dij{Q'^/Q'j), the factorization 
formula (I2.108p becomes 



dxdQ^ dxdQ^ 
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The scale Q'j is called the factorization scale, which we are free to choose. 
A typical choice is = cQ^, where c is between 0.5 and 2. Again, in all- 
orders calculation a physical cross-section does not depend on this choice, 
but in practice, as only the first few terms in the perturbative expansion 
are known, the prediction retains sensitivity to the adopted choice. 



2.6 Factorization for other processes 

Although 1 have here considered only the deeply inelastic scattering, the 
underlying physics is shared in variety of other processes involving hadrons 
in the initial state — the structure of the collinear singularities is the same. 
It follows that the parton densities should be independent of the actual hard 
process, universal. For example in the Drell-Yan production of dileptons in 
nucleon-nucleon colhsions, the leading logarithms originate from diagrams 
hke that in Fig. 12. 191 The formal proofs for factorization are highly technical 




Figure 2.19: Ladder- type diagram that gives rise to a leading logarithmic 
term. 

and mathematically demanding. Therefore, there are only few processes for 
which such all-order proofs actually exists [23j, but it is typically assumed 
that for hadronic interactions that are "hard enough" (involve a large in- 
variant scale), the pQCD-improved parton model is applicable. Ultimately, 
however, it is the comparison with experiments that is of essence. 
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2.7 Example of parton densities 



In their full glory, the factorization and DGLAP equations are exploited and 
tested in global QCD-analyses, to be described in ChapterEl In short, their 
purpose is to extract the x-dependence of the parton distributions from the 
experimental data. The Fig. 12.201 displays a typical outcome of such analy- 
ses, and Fig. [221] demonstrates how the Q^-dependence of the experimental 
F2 data becomes correctly reproduced by the DGLAP evolution. 
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Figure 2.21: Experimental data for proton structure function F2 from HI 
experiment [251 [26] compared with the CTEQ6M PDFs [S]. 
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Chapter 3 

Deeply inelastic scattering at 
NLO 



In this chapter I present the full calculation of the deeply inelastic scattering 
cross-section at next-to-leading order pQCD. Based on the discussion in 
Chapter O there will be three types of divergences: collinear, infrared, and 
ultraviolet. The most sophisticated and gauge-invariant way to regulate 
these is to use the dimensional regularization methods |27J and perform all 
calculations in = 4 — e space-time dimensions as in p8], or e.g. [29]. The 
divergences will appear as l/e"-poles but only those which correspond to 
the collinear divergences will remain. These are removed by absorbing them 
into the parton densities, after which the e ^ limit can be safely taken. 

The generic form of the hadronic tensor W^i, in Eq. ()2.6I) is independent 
of the space-time dimension, but the projections to the two independent 
structure functions Fi and F2 receive some e-dependence: 



F2 



X 2 - e 



MW^, (3.1) 



as may be directly verified, remembering that g^uQ^" 



3.1 Leading-order contribution 

For consistency, also the leading-order contribution must be computed in 
= 4 — e dimensions. The appropriate A^-dimensional extension of the 
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partonic tensor for the leading-order j*q q process in Eq. (j2.12p is 

AttMW^^, = f (2vr)iV-i2yo (2vr)^^(^^ {p + q - p') TV[/7>1 

= |^Tr[/7>1<5(e-x). 

The contractions with 5^1, and p'^p'^ give 

5^, Trb/'777'^] = -2(2 - e)^^ p^^p^ TVb/V|^7n = 0, (3.2) 
and we find 

- gf^uiMW^^,) = elz (1 - 1) 6{l - z) P^P^MWf,,) = 0, (3.3) 

_ 2 

where I have defined the variable z = x/^ = 
3.2 Gluon radiation 

The diagrams contributing to the real-gluon emission process 7*9 qg 
are shown in Fig. (13. ip . The squared, spin-independent matrix element 



p / 

Figure 3.1: The graphs for 0{as) real gluon radiation, 
contracted with —g^u reads 

= ^^glCF{^?rlH{l-^) (3.4) 

where gsifs^Y^^ is the strong coupling in = 4 — e dimensions and where 
the Mandelstam variables s, i, u are defined as 

s = {p + qf = {k + p'f 

u = {p — p')'^ = {q — k)'^ ■ 
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As we are eventually taking the e ^ limit, the last two terms in Eq. (I3.4p 
will not contribute and I will forget them from now on. The 2-particle phase 
space in N dimensions can be written as 



= (2'rf'"'^5+((p + <!-'-f) (3.5) 
In the center-of-mass frame of q and p, we may choose the momenta as 

p = (|p|,0,0,...,|p|) 

q = (<?0,0,0,...,-|p|) 

k = (|k|, . . . , |k| cos 6) 

p' = (|yO|,...,-|k|cos0). 

The squared matrix element does not depend on the momenta which were 
left unspecified above and we may perform the redundant angular integra- 
tions 

d^-^k = \k\^-'^d\k\nN-i (3.6) 

= \k\^-^sm^-^ed\k\d6dnN~2 



I 



d^N- 



where / diln is the surface area of an n-dimensional euchdean unit-sphere. 
The phase space (13. 5 p thus becomes 

dn(2) = i-JM!^(i|k||k|i^^d(cose) {1 - cos"^ ey"^^ 5{s - 2\k\VI). (3.7) 

The |k|-integration ehminates the remaining (5-function 

J d\k\\k\^-'5{s - 2\k\Vl) = ^(^^y , (3.8) 

and after introducing an angular variable y = ^{1 + cos 9), the 2-particle 
phase space becomes 
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In terms of y, Q^, and z = —q'^/2p ■ q, the Mandelstam invariants are 



^(1-^), i 



^ = -^y- (3.10) 



The square brackets from (J3.4p and y integral from the phase space (J3.9p 
leads to 



[y(i - y)] 



-e/2 



1 — z 1 — y 

\ 

l-y 1- z 



+ 2 



zy 



{l-z){l-y) 



which may be evaluated using a generic identity 



r(a + P) 



resulting in 

r^(l-e/2) 

r(i-e) 

Altogether, 



2-e 2z 
1-Z + 



1-zl 



+ 



1 - e/2 



2(l-e)(l-z) 



(3.11) 



(3.12) 



(3.13) 



Real 



(3.14) 



e2(l-e/2)Cp 



a. 



/4W/^r(l-e/2) 



e/2 



1 -Z + 



2^ V y r(i-e 

2z 1 \ l-e/2 



1-zl 



+ 



2(l-6)(l-z) 



where the designation "Real" reminds that this is a contribution from the 
tree-level real gluon emission. In this expression the collinear y ^ 1 diver- 
gences (gluon being emitted in the direction of incoming quark) are manifest 
as explicit 1/e-poles which now remain finite for e < 0. However, the ex- 
pression is also singular as z ^ 1 corresponding to the vanishing energy of 
the radiated gluon or gluon being colhnear with the outgoing quark. These 
singularities can be made exphcit by a distribution identity 



1 



log(l - z) 
1-z 



(3.15) 



(l-z)i+^ 

where the plus-distributions are defined as in Eq. (12.941) . Applying this 
identity to the lowest hne of Eq. ()3.14p gives 

" ^2 3 



1 



1 + z" 
1-z 

1 + 



2^(^ 



+ {l + z' 



7 



log(l - z) 
1 - z 



2 (1-z) 



logz + 3 - z + -5(1 - z), 



(3.16) 
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and we arrive at 



'^'^^ I Real « 



e«(i-^/2)c^i^^(^^j r(i _ 



5(1 -z) 



1 + z^ 
1 - z 



■<5(l-z) 



+ + 



log(l - 
1 - z 



-7 ^ logz + 3 - z + - (5(1 - z) ^ , 

2 1 - z + I- z ^ 2 ^ ^ ' ' 



(3.17) 



The contraction with P^^P^ is simpler and does not lead to singular be- 
haviour: 



Real 



1 



P^Pn-M7*.-..i;.rfn(2) (3.18) 



'^^ ' ^ 2-11 Az \xJ 



3.3 Virtual corrections 




Figure 3.2: Virtual corrections to deeply inelastic scattering. 

The virtual corrections to the Born cross-section Eq. ()3.3p stem from two 
sources: from loop correction to the photon-quark vertex and from the self- 
energy corrections to the external legs. These are shown in Fig. 13. 2[ In 
the Feynman gauge the 1-loop vertex correction (with massless external 
quarks) is simply a multiphcative factor [1 + T] to the tree-level Feynman 
rule —iej/j,, with 



r 



e/2 



r(i - 6/2) 
r(i - e) 



8 2 
- + 

e euv 



vr 



(3.19) 



where I have separately indicated the ultraviolet divergence euv > that 
occurs in the high end part of the loop momentum. The renormalization 
constant Zp for massless quarks is, in turn, 



Att \e euv 



(3.20) 
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The ultraviolet poles cancel, and the total virtual contribution is 



A**^^ Virtual 



/^^^^ Born 



\l + T\^Zl 



(3.21) 



el{l-e/2)5{l-z)^CF 



r(l - e/2) 



r(i-6) 



6 
e 



2 



2\ 



vr 
"3 



3.4 Total quark contribution 

When the virtual pieces and those from real gluon radiation are combined, 
the double poles evidently cancel, giving altogether 



{Mm,) 



A''^^ Virtual 



+ MWI) 



M'^^ Real 



r(i - 6/2) [ 2 

r(l-6) 1 e 



1 + ^' , 3^,^ 
— 0(1 — z) 



lofffl - z) 



2\ e/2 



+ + 



1 - Z 



1 



2 1-z 



i^log. + 3-.-(; + ^|*(l-.) 



The total quark contribution to the structure functions can now be obtained 
by folding them with the quark densities 



{MWf 



l^^'lBorn I Virtual 



^ Vi''' Real 



where the go(0 denotes the "bare", non-physical density which will even- 
tually go along the redefinition of the quark densities. When the various 
pieces are put together as instructed in Eq. (jS.ll) , 



1 



^2 



2 a. 



r(l - e/2) /47r/u' 



2X^/2 



where the coefficient function Cq{z) is defined as 



(1 + z^) 



log(l - z) 
1 - z 



1 + z^ 
' 1-z 



log 2; + 3 + 2z 



- 2(l-z)+ 

5 + 



(3.22) 



(3.23) 
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3.5 Initial state gluons 



The exact NLO calculation of the initial state gluon contributions proceeds 
much in a similar fashion as extraction of the quark contributions above. 
When averaging over the transverse gluon polarization states, one should 
remember that there are 2 - e such states instead of usual 2. The squared, 






Figure 3.3: The diagrams for 0{as) initial state gluon contribution. 



matrix element corresponding to the diagrams shown in Fig. 13. 3[ contracted 
with —g^,y reads 



(l-e/2)(| + ^)+2^^ 
t uj tu 



Supplying the phase-space element and doing the angular integrals, we end 
up with a following gluon contribution 



An" 

2 CKs / 4:TTfJ, 



27r V Q 



2\ 



2Tr 



r(l-e/2) [2.2/-, n2i 1 



1 - z 



The analogous results from the contraction with P^^P'^ are 



1 9-^qg\fj,iy 



(3.24) 



[z' + (1 - zf 



(3.25) 



pf'P'' MW. 



rG 



(iq-^^RQ 2 ' 



(3.26) 



From these the we can build the initial state gluon contribution to the 
structure functions 



X 



(3.27) 



2\ 
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where the coefficient function Cn(z) is defined as 



1 - z 



+ TR6z{l-z)-Pgg{z). (3.28) 



3.6 Total Fo and PDF-schemes 



We are now in a position to add up the quark and gluon contributions to 
the total F2. The remaining e-dependent terms come with a common factor 

,2\ 



2 /iTTfl' 



r(l-e/2) 

r(i - 6) 



jE + log(47r) + log 



(3.29) 



where I introduced a short-hand notation absorbing the Euler-Mascheroni 
constant and log(47r) to 1/e. Writing the total F2 explicitly, we have 

,2 ■ 



1 



X 



E 

q,q 



4lo 



1 



+ 2^62^0® I 



as 
2^ 

as 
'2tt 



Pqg + —C, 



27r 



Based on the discussion of Chapter [21 the collinear divergences, the 1/e 
pieces, should be absorbed into the redefinition of the quark and gluon 
densities which, as already noted, are not unique. In the dimensional reg- 
ularization framework, the general NLO definitions of the scale-dependent 
quark and gluon distributions at a certain factorization scale can be 
written as 



Ax,Q)) 



9{x,Q)) 



go® < 1 



50 O < -7;- 



50 ® < 1 




a 



s /■scheme 



q,q 




p I !■ scheme I 

I 

( 

1 

p I scheme I 

^<;5 + 2^S J 



(3.30) 



P 



99 



scheme 

27r^^ 



where /|^^'^™'^, h^^^^'^'^ are arbitrary (finite) functions. It should be under- 
stood that these terms are just the first ones of a whole series which formally 
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exponentiate and give rise to the DGLAP evolution as discussed in Chap- 
ter [2l With these definitions the expression for the NLO structure function 
F2 becomes 



1 



a 

2tt L 



scheme 



+2Y^elgiQ})^\ log ( g ) P.. + ^ ^ - 



2tt 



^scheme 
'9 



(3.31) 



where 



log(l - z) 
1 - z 



1 + z^ 
1-z 



loff 2; + 3 + 2z 



Cgiz) = Pggiz) lo: 



1-z 



1 



- 2 (1 - z)+ 

2 + Yr(^-^' 



+ TR6zil-Z)-Pggiz). 



(3.32) 



(3.33) 



Obviously, there are arbitrarily may ways to define the scheme — the two 
most common ones are 

• MS-scheme: 

In this scheme only the colHnear divergence and the regularization- 
framework-originating — 7£; + log(47r) are absorbed into the definition, 

i.e. one chooses /«=heme ^ ^scheme ^ ^^scheme ^ ^scheme ^ q. 



DIS-scheme: 



Cr, 



This scheme is defined by choosing /^^^^heme ^ ^scheme _ 

which maintains the simple naive parton model form of the cross- 
section at = Qj. How to define /i«=heme ^^^j /^scheme DIS-scheme 
is, however, more or less a matter of convention. Usual choice ^ is 
dictated by the momentum conservation 



MS 



q,q 



X 



DIS 



which requires M^heme ^ -2nfCg and M^heme = _Cg. 
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For completeness, I record here also the expressions for Fi: 



q,q 



as ^ 



^r«2z(i - z) 



(3.34) 



3.7 Numerical estimate 



The numerical calculation of deeply inelastic structure functions and cross- 
sections at NLO is relatively simple with only single integrals to be numeri- 
cally evaluated. In order to get in touch with the expected size of the NLO 
corrections, I have computed the proton F2 both in the leading and in the 
next-to-leading order (MS-scheme) with the same set of parton densities 
(CTEQ6L1 [5T]). The Fig. [Ml presents the results at two typical scales, at 
Q2 = lOGeV^ and at = 100 GeV^. In most x- values the expected NLO 




Figure 3.4: The left panel shows the behaviour of absolute proton F2 at = 
lOGeV^ and = 100 GeV^ (with offset +0.2) calculated with CTEQ6L1 
leading order parton densities in the leading-order and in the next-to-leading 
order. The panels on the right display the NLO contributions from the quark 
and gluon channels separately normalized by the leading-order F2. 

corrections are of the order of few percents and only the quark contribution 
shows a growing behaviour at very large x. Although the leading-order term 
dominates and F2 as such is not very sensitive e.g. to the gluon content of 
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the nucleoqj, such terms cannot, however, be neglected in a precision anal- 
ysis. An example of a quantity that is directly much more sensitive to the 
gluons is the longitudinal structure function Fl = F2 — 2xFi which vanishes 
at leading order, as can be seen from Eq. 13.341 

Also the NNLO coefficient functions for the deeply inelastic structure func- 
tions are nowadays known [32], and an analogous estimate as shown here 
reveals the expected magnitude of the NNLO vs. NLO corrections similar 
to the NLO vs. LO shown here. In other words, the NNLO terms in the 
coefficient functions are still important. 



^Gluons, however, are the driving force in the DGLAP evolution of small-a; quarks and 
thus their effect on F2 is significant but indirect. 
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Chapter 4 

Drell-Yan NLO cross-section 



The Drell-Yan dilepton production in nucleon-nucleon collisions is another 
process which can be employed in probing the parton distributions. Here, 
1 will derive the double-differential NLO cross-section (Pa /dM'^dyR, where 
and hr denote the invariant mass and the rapidity of the lepton pair. A 
detailed computation of this particular cross-section has proven to be sur- 
prisingly difficult to find from the literature: Up to my knowledge, such can 
only be found in [33] which employs a massive gluon scheme to regularize 
the colhnear and infrared singularities. However, in order to rigorously em- 
ploy the NLO definitions of the scale dependent parton densities discussed 
in the context of deeply inelastic scattering in Chapter 2, the dimensional 
regularization methods are to be used. The results obtained in this section 
coincide with those given in [34] . 

4.1 Leading-order calculation 




Figure 4.1: Leading-order parton diagram for Drell-Yan dilepton produc- 
tion. 

In the leading order, the lepton pair production proceeds by annihilation of 
quark and antiquark as shown in Fig. 14.11 The unpolarized, color- averaged 
square of the matrix element corresponding to this diagram can be com- 
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pactly written as 

12q 



= Tra^'^'^' (4-1) 



where, in N dimensions, 



L'^- ^ Tr[|^i7'^^27l (4.2) 



^ Tr 



When the leptonic tensor L^^'^ is integrated over the 2-particle phase space, 
the dependence on momenta ki and ^2 washes out and the result only 
depends on the virtual photon momentum q. Moreover, since q^L^^ = 
quL^^^ = 0, the Lorentz structure is necessarily 

j dn(2)L^'^ = (^gi^^ - L(g2), (4.3) 

where 

Similarly, q^'H^^^ = q^H^^ = and therefore 

j dU^^^L'^^ = [g^-'Hl,] L{q^). (4.4) 

Thus, the cross-section splits into independently calculable leptonic and 
hadronic pieces 

al = f dn(2)L'^-i/0 (4.5) 

° 2s 12g4 J ^ ^ 

1 e^el , „„ „^ , 1 



2s 12g4 

simplifying the calculation. Prom ()4.2p , we have 

g^^L^'' = -2(2-6)g2 (4.6) 
g^'^H^, = -2{2-e)s. 

Since the leptonic quantity J dll^^) {g^^L^^'^) will be common also for 
the higher order diagrams, it is unnecessary to drag its exact A^-dimensional 
expression throughout the calculation — what matters is the hadronic part 
g^'^Hl,. In the iV ^ 4 hmit / dlP^ = ^, and (jH]) gives 

*2 = ^^?. (4.7) 
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The invariant mass M^, and the rapidity yu of the produced lepton pair are 
defined as 

M2 = {ki+k2f=q^ (4.8) 

At leading order, fci + /c2 = pi +P2, and in the center-of-mass frame of the 
colhding nucleons 

Pi = ^6(1,0,0,1), p2 = ^6(1,0,0,-1), (4.10) 
it follows that 

1 f 

M^ = iii2S, yK = -log|i. (4.11) 

Thus, the double-differential partonic cross-section in these kinematic vari- 
ables can be written as 



dM^dyn 



al5{M'' - 66s)<5 [vR - ^ log I) (4-12) 



Ana, 



2 

em 2 



9sM2 



fig 5(1 - 2;)(5 [vR - ^log^) • 



where z = M"^ /s. To obtain the corresponding hadronic cross-section we 
integrate over the parton densities of the incoming nucleons and sum over 
all flavors 

Performing the integrals we find 



dW^dyR 9s 



where 

xi = A/re^«, X2 = ^/^e"^«, = — • (4-13) 

4.2 Gluon radiation 

The QCD corrections to the Born-level cross-section from emission of one 
additional unobserved gluon are shown in Fig. 14.21 The corresponding par- 
tonic cross-section for such process can be written as 

al = 11 1 e'^el f dn^^)^^-//^^, (4.14) 



^ 2s 2_><_2 X 3_x_3 ./ 

spin color 
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Figure 4.2: Gluon radiation in Drell-Yan dilepton production to NLO 
pQCD. 



where L^^^ is the same leptonic tensor as in the previous section, and Hj^^, is 
obtained from the spin- and color-summed square of diagrams in Fig. 14.21 
In order to separate the leptonic and hadronic pieces as was done in the 
leading order, we use an identity 



1 



2-K 



(27r)^-i2£:, 



(fci+fc2)2=M2 



to rewrite the 3-particle phase space as 



(27r)^-i2i?fc, (27r)^-i2ii;fc, (27r)^-i2i?fc3 
{2^f5^^\pi+P2-ki-k2-h) 



(4.15) 



(4.16) 



d^~^ko 



{27T)^-^2Ek, {2t,)^-^2Eu, 



{27:f6'^^\q-ki-k2 



d^-^k. 



d^- 



(27r)^-i2^fc3 (2^)^-12^, 



■(27r)^<5(^)(pi+P2-fe3-g) 



dM"^ 



2tt 



If- 



By this trick, 

daj^ 



1 1 



dM2 



72sM4 27r - 1 



1 



1 



9sM2 67r 2(2 - N) 



dli^^^^'^H^ 



(4.17) 



where I already took the e ^ limit of the leptonic part f dIi^l\^uL^^'" 
In the present case, we define the Mandelstam variables 

M2 



u 



{pi - hf 

iP2 - ksf 



M2 

-Mia 

z 



il-z){l-y) 



z)y, 
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where I have introduced the angular variable y = ^(1 + cosO) G [0, 1], with 
9 referring to the angle between incoming quark and emitted gluon in the 
center-of-mass frame of the incoming quark and antiquark 

pi = (Vi/2,0,0,... , Vi/2) 

P2 = (\/I/2,0,0,...,-\/I/2) 

h = (|k3|,...,|k3|cos0); \ks\ = {s-M^)/{2Vl) 

q = ((?°,...,-|k3|cos0); = (s + M2)/(2VI) 

and z = ]VP/s G [r, 1]. The rapidity of the produced lepton pair in this 
frame is 



■ log ■ 
2 + 



(4.18) 



and owing to the additivity under Lorentz-boosts, the rapidity ur in the 
nucleon-nucleon center-of-mass frame is 



VR 



- log — + = - log 
2 ^2 2 



Ci 1 - y(i 



_^2Z + y{l-z)_ 

Thus, the wanted doubly differential cross-section is obtained by 



1 



VR 



■log 



ii l-y{l-z) 
i2Z + y{l-z) 



(4.19) 



(4.20) 



Having now fixed the kinematics, we can write down the squared matrix 
element, contracted with the metric tensor g^'^, 



gf^^Hf^, = -gH2CF{l-e/2) 



(2 



u t . 
^ + - 1 + 4- 

t u , 



til 



(4.21) 



In terms of the variables z and y this becomes 



-gH2CF{2-e) 



(1 - 6/2) 



1 



V 



+ 



+ 



2z 



A similar calculation that led to Eq. ()3.9I) gives the 2-particle phase space 



dU 



(2) 
H 



87rr(l 



47r 
M2 



e/2 



dy - y)] 



-e/2 







(4.22) 

At this moment, the various poles that occur in the kinematic corners z ^ 1, 
y ^ 1, y ^ should be made exphcit by using the following distribution 
identities 



1 



1 



(l-z)i+^ 
1 



5(1 -z) + 



1 



yl+e 



--5(z) + - 

6 Z+ 



z) + 
log z 



logfl 
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After some algebra, we find the following stack of terms 

(1 - 



(1 - 6/2) 



1-1/ ^ y 



5{l-z) [5{l - y) + 5{y)] + - 



+ 



2z 



e2 e 



+ S{l-z) 



1 1 

+ 



e L(i-y)+ y+. 
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All terms in the third hne above are irrelevant as they vanish under inte- 
gration: for example 
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We also notice the following analytic result 
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-C(2)= 



■y<^(y«-^log|) 5(1-^), 



where C(2) refers to the Riemann zeta-function. Performing the remaining S- 
function-restricted y-integrals, we reach the final form of the partonic cross- 
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section: 
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where I have defined 
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(4.25) 



4.3 Virtual corrections 



The virtual 1-loop diagrams are essentially same as for deeply inelastic scat- 
tering, Eqs. ()3.19p and ()3.20p , but the vertex contribution must be analyti- 
cally continued to the time-hke region = —q"^ — < 0. The virtual 
contributions may be written as 
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Adding the virtual and real gluon emission contributions, the infrared 
poles cancel, giving 



dyndM^ dyndM^ | Victual 9sM2 
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To turn the partonic cross-section above to a hadronic one, we integrate 
over the parton densities 
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'2^9 
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^o(6,6), 



where 



The 0-subscript above is to remind us that these are still "bare" parton 
densities, anticipating their replacement by the scale-dependent ones. Per- 
forming the integrals constrained by the delta-functions, we obtain 
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Figure 4.3: Gluon-quark subprocess in Drell-Yan dilepton production. 



4.4 Initial state gluons 



The second type of 0{as) process that contributes to the Drell-Yan cross- 
section is the one with a gluon in the initial state. The two diagrams for 
such process are shown in Fig. 14.31 The corresponding partonic cross-section 
can be written as 
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and it follows that 
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The hadronic part RfS is the color- and spin-summed square of the diagrams 
shown in Fig. 14.31 
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In terms of variables y and z, defined as in the previous section, this reads 
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The phase-space is 
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In this case only collinear y ^ singularities are present, which can be 
turned into explicit 1/e-poles by distributions: 



b(l-y)]-/2(l-z)i-V/2<{ (1-6/2) 
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Putting all factors together, we have 
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the hadronic cross-section becomes 



where 
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The contribution from the mirror process is obtained similarly. Defining 

K^\iu^2) [g'i\i,) (^f (6) , (4.36) 
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the result is 



(4.37) 



2z{l - z) + (1 - z)\l - yo) + 



z' + {l-; 



4.5 Finite result 



The results derived above still contain all collinear 1/e-poles. The uni- 
versality of the parton distributions requires that it must be possible to 
remove these singularities by the same definition as has been apphed in 
deeply inelastic scattering. To first order in strong couphng we can invert 
the definitions (j3.30p to write 



Ho{xi,X2) = H{xi,X2,Q}) (4.38) 




and 

Kf^'''{xi,X2) = K''''''"'{xuX2,Q}) + 0{as). (4.39) 
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Inserting these expressions to ([09]) , all 1/e-poles in (109l) , ([05|) and ([QTll 
cancel giving the final, finite results: 
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It is not very difficult to integrate these expressions over y/j to recover the 
differential da/dJVP cross-section given e.g in [48j. 



4.6 Numerical implementation 

The presence of double integrals makes the numerical calculation of Drell- 
Yan cross-section somewhat more challenging than the deeply inelastic scat- 
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tering. Especially, one should pay attention how to evaluate integrals in- 
volving a product of plus-distributions. First, by change of variables 




Figure 4.4: The integration regions for hr > 0. The constant yo-lines and 
z-hyperbola are indicated. 



Jx^ ?1 Jx2 ?2 ^ 1 - 2 V ^ « / J Ja+b+c z 

where the division of the rectangular integration region in (^i, ^2)-plane to 
A, B and C parts are indicated in Fig. 14.41 for ur > 0. Among other terms, 
the NLO contribution to quark-antiquark process involves a term 

1 + ^2 



(l-z)+(l-yo)+' 
According to the definition of the plus-distributions, 

dyoT^ ^T, ^ = / —T, ^ / dyo- 
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where 

lz + y{l-z) T 
\] zl - y[l - z) z^i 
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Applying the definition of plus-distributions again to the remaining z-integral, 
we have 



Jo (T^io 
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where the subtraction term evidently vanishes, leaving 
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The regions B and C avoid both the z = 1 and the y = 1 singularities, 
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Since H{^i,(2) = for ^1,^2 > 0, the j/o-integrals for regions B and C may 
be extended to range from to 1, and all integrals can be grouped neatly 
together: 
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Although we considered here the special case yu > 0, the result above is 
vahd also for yR < 0. The treatment of other integrals involving plus- 
distributions is a straightforward extension to what was presented above 
and for completeness, 1 record here the M5-scheme cross-sections in a form 
which no longer explicitly involves plus-distributions and can thus easily be 
turned into a computer code: 
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4.7 Numerical estimate 
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In order to get a feeling about the size of the NLO corrections to the lead- 
ing order cross-section, Fig. 14.51 shows a numerical result in the MS-scheme 
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Figure 4.5: An example of the rapidity distributions in Drell-Yan process 
computed in the MS-scheme. The left panel shows the behaviour of the 
absolute cross-section for ^/s = 38.76 GeV and M = 8 GeV using CTEQ6L1 
leading-order parton densities and indicating the leading order and next- 
to-leading order contributions separately. In the right panel, the relative 
contributions normahzed to the leading order cross-section ctq of quark- 
antiquark and (anti)quark-gluon processes are shown. 

with the factorization scale choice Q'j = M^, for a kinematical configu- 
ration ^/s = 38.76 GeV (corresponding to a fixed target experiment with 
800 GeV proton beam) and M = 8 GeV using the CTEQ6L1 leading-order 
parton distribution functions. Unlike in the deeply inelastic scattering, the 
contribution of the NLO terms is rather large — always at least 50% and 
increasing when going away from the midrapidity. It is interesting to notice 
that the quark-gluon contribution remains always negative, partly cancelling 
the large positive contribution of quark-quark subprocesses. 

The Drell-Yan rapidity distribution, computed here to NLO, is nowadays 
known still one power higher in as (NNLO) [34]. The size of the NNLO 
corrections relative to the NLO are not, however, as large as the NLO relative 
to LO and the perturbative expansion seems to stabihze. 
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Chapter 5 

Inclusive hadron production 



For a deeply inelastic scattering event the experimental signature typically 
consists of the scattered lepton and a narrow shower, a jet of hadrons, orig- 
inating from the struck parton that balances the transverse momentum of 
the scattered lepton. The hadronization — how the partons transform 
themselves into a cascade of hadrons — is a non-perturbative process and 
beyond the reach of pQCD-tools. As long as we are not concerned about the 
structure of the jet, we can simply ignore such process. This kind of final 
state is said to be fully inclusive. However, as the jets consist of hadrons, 
with a suitable detector they can be identified and their momentum mea- 
sured. In this Chapter 1 shortly describe how such indentified hadron pro- 
duction cross-sections are calculated in pQCD, especially in hadron-hadron 
collisions. 

5.1 Fragmentation functions 

Let us return to the leading-order deeply inelastic scattering, where a high- 
photon knocks a quark q to an escape-course from the nucleon triggering 
off a jet as shown schematically in Fig. 15.11 




Figure 5.1: Hadron production in deeply inelastic scattering. 
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The number density of hadrons h carrying a fraction z of the jet energy is 



described by a fragmentation function Dg{z), where the q indicates the 
parton that initiated the fragmentation. Consequently, the leading-order 
cross-section for single particle (plus anything else) production in the deeply 
inelastic scattering (SIDIS) is 

- '"^E^im^^iz), (5.1) 



dxdzdQ"^ dxdQ"^ 



q 



where the energy-fraction z may be expressed in an invariant form as 

P ' P -E'hadron 



(5.2) 



where the latter equality refers to the target rest frame. Beyond leading 
order, however, collinear divergences due to e.g colhnear gluon radiation 
from the outgoing quark emerge. Here, I would emphasize that these diver- 
gences remain only because the final state is not inclusive enough: If it was 
not for the fragmentation functions — if we would not care what the quark 
will eventually become of — the divergences above would exactly cancel 
against the loop diagrams. In the axial gauge, the dominant 0{as) loga- 
rithm originates from the quark-sphtting diagram shown in Fig. 15. 2[ and 





Figure 5.2: Dominant gluon radiation graphs for deeply inelastic hadron 
production. 



the contribution to the cross-section can be expressed as 

where again denotes the convolution integral and the splitting functions 
Pqq and Pgq are the same as earher. These leading logarithmic terms can 
be resummed essentially in the same way as was done in Chapter [2] for 
the initial state radiation. Then, by absorbing these logarithms into the 
definition of scale-dependent fragmentation functions D'^{z,Q'^), they are 



dxdQ^ ^ ^ 



2tt 
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seen to follow the same DGLAP equations 
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ag2 


27r L 



Pm, ® K iQ') + Pm ® D^giQ^)] (5.3) 



Pgg®D^^{Q'')+Pg,^(^Dl{Q 



as the parton distribution functions do. At the leading logarithmic ap- 
proximation, the sphtting kernels are exactly the same as derived in the 
initial state parton branching, but at higher orders they become different 
being still related by a proper analytic continuation [35]. As was argued 
in Chapter O the collinear divergences between the incoming and outgoing 
partons do not interfere. That is, they independently factorize, and the first 
pQCD-improved version of Eq. (j5.1l) would be the one in which the parton 
distributions and fragmentation functions are scale-dependent. 



dxdzdQ^ dxdQ"^ 



Y,elf,{x,Q^)D^iz,Q^ 



(5.4) 



For an exact NLO calculation, one should consider the lepton-parton pro- 
cesses 

i{k)+p\^4 + e{k') + x, 

and adopt the dimensional regularization methods to calculate the differen- 
tial cross-sections 
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-*P'2- 



dxdzdQ^ 



2p\ • q 



P-q' 



The diagrams for such calculation are the same as those displayed in Chapter 
m but the phase-space integrals should be constrained by an additional delta 
function 5{z — {P ■ P2)/{P ■ q)). I will not get into details here |M1, but at 
the end, the partonic cross- sections will retain several 1/e-poles from the 
colhnear divergences. Multiplying these expressions by the "bare" parton 
densities and fragmentation functions one can construct the hadronic cross- 
section 



dxdzdQ"^ 



dxdzdQ'^ 



^0 



(5.5) 



where p\ = ^iP and = P' 1^2- Defining the scale-dependent fragmenta- 
tion functions similarly to the parton distributions in Eq. ()3.30p , and writing 
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the bare quantities fio{x) and Dq{x) in terms of the scale-dependent ones, 
one can take the e ^ hmit 




1 



dxdzdQ"^ 



dxdzdQ^ 



where da depends on Q"^, Qf^^^, and Ql^^,. This is a cross-section guaranteed 
to be free from colUnear divergences. 

The fragmentation functions can be determined by analyzing experimental 
data by a similar procedure which is used to extract the parton densities (to 
be described later) |33 [38l [Ml [40] . The cleanest enviroment to extract the 
fragmentation functions are the e+e~-induced processes due to less crowded 
final state (compared to the colhsions involving hadrons), and as they are 
not complicated by the parton distribution functions. However, the e+e~- 
data alone cannot constrain all components of the fragmentation functions 
well, and therefore some latest analyses like [39j are complemented by SIDIS 
measurements just introduced, and also by data from hadron-hadron col- 
hsions [Ml [Ml [iQ] to be described shortly. The resulting fragmentation 
functions obviously somewhat depend on the particular set of parton den- 
sities used in the analysis, and the most complete analysis would combine 
both to a single analysis. 

5.2 Single-inclusive hadron production in hadronic 
collisions 

The fragmentation function formahsm can also be apphed to inclusive pro- 
duction of large transverse momentum hadrons in hadron-hadron collisions: 



where Ki and K2 denote the incoming hadron momenta, is the momen- 
tum of the observed hadron, and X is, as usual, anything. The calculation of 
this cross-section begins by computing the invariant partonic cross-sections 



At NLO [H], these consist of three pieces: purely tree-level 2^2 and 2^3 
diagrams, and 2^2 diagrams decorated with loops. As earher, there will 
be various divergences: colhnear, infrared, and ultraviolet. The singularities 
appearing as 1/e^-poles cancel between the real and virtual contributions, 
and when the ultraviolet divergences from the loop integrals are subtracted 
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(5.6) 
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according to the adopted renormalization prescription, only the collinear 
1/e-poles remain. Multiplying these quantities by the (bare) parton densities 
and fragmentation functions and integrating over the available phase-space, 
the invariant hadronic cross-section becomes 
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where the additional factor 



originates from d^K^/E^ = xld^p-^/p^. 
Writing the bare quantities fio{x) and Dq{x) in terms of the scale-dependent 
ones, the remaining colhnear divergences again cancel and a finite result is 
obtained in the e ^ limit. 
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The cross-sections are often reported specifying the center-of-mass energy 
of the hadron-hadron collision, the transverse momentum pr, and the 
rapidity y^, of the observed hadron. In terms of these variables, the integra- 
tion limits in the expression above are 
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x\pTe 



-VR 



2pT cosh yR 



prpe vr ' " x\^fs — ptg^i^ \Jx\X2s 

Similarly to the Drell-Yan dilepton production, the leading-order approxi- 
mation to the inclusive hadron production turns out to undershoot the ex- 
perimental data [I2] by a typical factor of ~ 2, depending on the kinematical 
variables and scale choices. The larger cross-sections at NLO improve such 
situation nicely although at low y/s < 60GeV the theory still seems 
to undershoot the data. For more recent work discussing RHIC data for 
constraining the fragmentation functions, see e.g. [38] |39] . 
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Chapter 6 

Solving the DGLAP 
equations 



Global QCD-analyses, to be described later in Chapter [71 require an efficient 
way of solving the DGLAP evolution equations. Being integro-differential 
equations, there is not much that can be done purely analytically but nu- 
merical methods are to be used. Several methods to accomphsh this has 
been developed — for an elementary account for couple of treatments at 
leading order, see [H]. At leading order, the DGLAP equations are still 
fairly easy to solve but the technical difficulties significantly increase when 
going to higher orders (NLO & NNLO). The method I describe in this sec- 
tion is based on |45] and it has been employed in the pubhcations [llH [iVj 
of this thesis. For description of further methods and available codes see 
e.g. [MliTj. 



6.1 Decomposition of the DGLAP equations 

The full set of evolution equations to be solved can be written as 



Q 
Q 
Q 



2 dqi 


_ as 






2 9qi 


_ Ols 






2 dg 


_ as 




2tt 



Yl ^I'lk ^Qk + Y, Pi^Qk ^ Ik + Pqg(^9 
. k k 

Y ^Qk + Y Pq^qk (^Qk + Pqg^g 

. k k 



(6.1) 



where the arguments of parton densities and the strong couphng are not 
displayed. A useful decomposition [I9l HH] of the sphtting functions Pqq 
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and Pqq is to separate the flavor-preserving "valence" and possibly-flavor- 
changing "sea" parts as 
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At leading order only P^^ in (I6.2p is non-zero, but at NLO they all are 
non-trivial, but respect the following relations: 
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which are reflections from the charge-conjugation invariance and the SU(3) 
flavor symmetry, but can also be easily understood on the basis of the Feyn- 
man diagramqj. By defining 
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where n/ is the number of active flavors, and 



(6.4) 



Qi ± Qi 



± 



(6.5) 



the set of equations ()6.ip can written as 
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(6.6) 

(6.7) 
(6.8) 



The densities q~ and qf — evolve independently, whereas and 

g are coupled. The strategy to solve the evolution of individual flavors qi, 
is to substitute q~^ derived from ()6.6p to result of (16. Sp and use (16. Sp . A 
good reference containing the expressions for all splitting functions needed 
to solve (lM])-([63]) is [49j. 

^At NNLO, however, P^q / P,^ 
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6.2 The Taylor expansion 



To keep the subsequent discussion as transparent as possible, let us con- 
sider the simplest evolution equation, namely that for the valence quarks 

q^{x,Q'^)=q-{x,Q^) (with P = P_), 



(6.9) 



27r 



■P ® qv, 



with a given initial condition qvix.Q"^). To make the Q^-evolution appear 
as hnear as possible, it is useful to define a new evolution variable 



(6.10) 



where /3o = ^Cg — ^^RUf appears in the QCD renormahzation group 
equation 



m — -. 1- Pi 

47r 
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(6.11) 



Trading the Q^-derivative with t-derivative, we have 



^2 d ^ asiQ^ 
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2/3o 27r 



+ 0(ai 



(6.12) 



With this change of evolution variable, the Eq. (16. 9 p reads 



1 + 



(3i as{Q^) 
2/3o 27r 



dt 



qv{x,t) = P®q.u{t). 



(6.13) 



To the NLO accuracy, the sphtting function P is of the form 



P(z) = pW(z) + ^P«(z), 
27r 



(6.14) 



and we may write Eq. ()6.13p as 



d 

■^qv{x,t) = O.0q^{t), 



(6.15) 



where 



= + 



27r 



p(l) _ p(0) 

2/3o 



(6.16) 
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The very crux of the matter here is to expand qy{t) as a Taylor series around 
the initial scale to = t((5o) = 

qv{x,t) = }_^- , (6.17) 

*;=0 

where d^qv{x,Qi)/dt'^ are multiple derivatives 

dqvit) 



dt 



Qvit) 



d^qJt) d^n , , agJt) dVL d'^qJt) ^ d^qJt) 



By using the lower-order derivatives in the expression for the higher deriva- 
tives, the nth one we can be written as 

^^ = MW®g„(0, (6.18) 

where each M^'^^ can be iteratively computed from previous ones 

= 1 

M(2) = 1^(i) + q(o)0mW 

M(3) = ^ 2^(1) m(^) + m(2) 

M^'') = j^(^) + 3J^(2) ® m(^) + 3^^(^) m(2) + j)(o) ® m(^) 



where = and O^*^) = = 0)/cii*^ for k>l. In general, 

jy-W = ^ ('^"^')j^W®m('=-i-") (6.19) 

«=o V / 



j^(0) ^ p(0)^ «.(Qo) (p(i)_i^p(o)) (6.20) 
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n(n) ^ (_^)""-^(Qo) (p(l)_i^p(0)) (6.21) 
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where [f) = f^n^^Lf^y is the usual binomial coefficient. Thus, the Taylor series 
(I6.17P becomes 



qv{x,t) 



.k=0 



(g) qv{x,0). 



(6.22) 



The crucial feature to be noticed is that the for each x, the evolution func- 
tions are independent of the parton density q^. Also, the magnitude 
of t in the physically conceivable domain is ~ 0.1, and one can expect the 
series to converge with a reasonable number of terms in the expansion. 



Since the parton densities tend to generally diverge as x~^, {6 > 0) towards 
X ^ 0, it is numerically more stable to damp such bad behaviour by writ- 
ing the evolution not for the absolute parton density q{x,Q'^) itself, but 
rather for the momentum distribution q{x,Q'^) = xq{x,Q'^). Multiplying 
the expansion (16.221) by x, 



qvix,t) 



^ A;! 



.A;=0 



(g)gt,(x,0), 



(6.23) 



where the "hatted" convolution should be understood as 



/l®/2 



fN{CN)q{ 



6.3 Integration 

In order to actually calculate the very formal Taylor expansion written down 
in the above section, one needs to evaluate a series of "nested" integrals, each 
one of the general form 

jjmonji (6.24) 

where F is a result from a similar integral. To accomphsh this task it is 
useful to break the x-interval [x, 1] into N smaller sub-intervals by a dis- 
crete grid (xq = x,xi,X2, . . . ,XAr_i,XAr = 1), and approximate the function 
-F(x) in each interval by a simpler one, like a polynomial, for which the 
integrals against the sphtting functions can be analytically evaluated. In 
other words, we write 

m 

F{x) ^^af\ii{x), Vx G [xfc,Xfe+i], (6.25) 
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where o'^'^ are coefficients which naturally depend on F{x). For example, 



if we employ a 3rd order polynomial, gi{x) 
aj^'^ can be taken to satisfy 



m 



4, the coefficients 
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(6.26) 



V "4 / 



that is, the coefficients of the polynomial are chosen to match the F{x) at 
four points around Xk as illustrated in Fig. ()6.ip . By inverting (16.261) , 



■,F{x] 



Xk-l i-l+i Xty, Xk+3 



Figure 6.1: Illustration of a polynomial fit to the function F{x) used in the 
interval x £ [xk,Xk+i]- The black dots denote the matching points. 
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(6.27) 



The various terms in the splitting functions P can be grouped to the fol- 
lowing categories 



P{z) 



A{z) 



(1-z). 



+ B{z) + C5{l 



(6.28) 
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and inserting such expression to (I6.24p , we find 





dz ' 




= Xi — 







+ 



Z — Xi 2^ 



C + A{l)logil - Xi) F{xi). 



Decomposing the integrals above as = Ylk=i Ix^^^ ' using the ap- 
proximation (I6.25P , we find 



m N~l m 

S E«fH72+p2) (6.29) 

e=l k=i+l £=1 

[C + A{1) log(l - X,) - A{l)a'] F{xi), 

'+1 dz A{xi/z)ge{z) - A{l)gi{xi) 



where 
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(6.30) 
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Xi 



1 dz 1 



Substituting here the coefficients from Bq. ()6.27p , the integral 
can be written as 

\ ' Jxi S ._n 



fe=0 

where the entries of the matrix Pik can be computed from 

4 4 



r=l i=l 

4 4 N-l 

r=l e=l n=i+l 

+ [c + A{l)log{l-Xi)-A{l)a 



k=i + T-2 



(6.32) 



k=n+r—2 



k=i 



In this way, the sphtting functions P and consequently also the functions 
fi^") and M(") derived from those, become matrices and the multi-dimensional 
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integrals reduces to mundane matrix multiplication 



oo N 



,k=On=0 



kV 



(k) 
in 



(6.33) 



where the matrices M do not depend on the form of the parton density, but 
only on the x-grid, and the initial scale Qq. That is, they can be computed 
once and for all. The hard manual work is to analytically evaluate the inte- 
grals in Eq. (16.301) for all splitting functions. At leading order the expressions 
are still reasonable to compute even by hand, but already at NLO-level the 
expressions become long — involving special functions hke polylogarithms 
and Riemann zeta functions — and use of a symbolic computer program 
like Mathematica is, in practice, mandatory. 



6.4 Numerical test of the method 

Using the NLO sphtting functions given in [49] (and also the leading-order 
ones discussed in Chapter [2]), 1 have constructed a Fortran code for calcu- 
lating the evolution with the method described above. In order to verify 
the accuracy of the method, 1 have tested it against the benchmark values 
of parton distribution functions given in Ref. f50]. This reference contains 
carefully cross-checked results for evolution of a given initial parametrization 
of partons from 2 GeV^ up to 10000 GeV^ with an unambiguously defined 
as- For doing the comparison, 1 have constructed an x-grid from x = 10~^ 
to X = 1 with 100 logarithmic intervals in a range 10^^ . . . 10~^ and 100 
hnear intervals in a range 10~^ ... 1. I have truncated the Taylor expansion 
to include the first 9 terms. The results for specific combinations of partons 
are displayed in Fig. 16.21 as a relative difference to the benchmark partons 
in percents. The evidently excellent agreement with the benchmark partons 
proves the apphed integration method accurate as well as that already the 
9th order Taylor expansion seems to convergence nicely. Only at very large 
X, where the parton densities are numerically very small — typically falhng 
off like (1 — a > 3 — a third-order polynomial does not optimally fit the 
input densities and a higher order expansion would be needed. However, for 
the purposes of our global PDF analyses the very large-x region is rather 
irrelevant for all other parton types except maybe the valence quarks which 
are numerically larger. 
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Figure 6.2: Comparison to the benchmark NLO partons [50j for valence 
quarks, sea quarks, and gluons. 
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Chapter 7 

About global QCD analyses 



Within the pQCD-improved parton model, the hadronic cross-sections for 
hard scattering processes can be calculated through the factorization theo- 
rem folding the universal, scale-dependent PDFs fi{x,Q'j) with the pertur- 
batively calculable pieces da, formally 

The physical content of this expression has been discussed in the preceding 
chapters. Thus, experimental measurements provide information about the 
PDFs fiiQ'j) as well as about the underlying parton dynamics da. This is, 
in short, the central idea behind the global QCD-analyses of PDFs. Here 
the word "global" is related to the universahty-hypothesis of the PDFs i.e. 
their process-independence: As much experimental data as possible should 
be considered simultaneously to find whether this is really true — can they 
all be described by same set of PDFs. If not, it may be a sign of factorization 
breakdown or perhaps discovering new physics beyond the Standard Model. 
However, due to the enormous complexity of the present-day accelerator- 
based experiments, one should also be cautious of not being misled by e.g. 
fluctuations in the data that might not respect any textbook statistics. It 
should be emphasized that the global analyses do not only constrain the 
PDFs, but also various fundamental parameters hke the strong couphng Os, 
the heavy quark masses, and even the elements of the CKM-matrix. 

The modern global analyses employing data from several free proton ex- 
periments was ushered in by works of Morfin and Tung triggering 
an enormous effort which today demonstrate a huge success with continu- 
ously increasing amount of data accommodated in the analysis. The leading 

^ Sadly, Wu-Ki Tung passed away during the writing of this thesis in March 2009. 
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groups in this domain are nowadays the CTEQ [52], the MSTW [53], and 
the NNPDF [54] collaborations, but various other parties like the HERA- 
PDF consortium (see e.g. [55]) which often focus only on a more restricted 
data input, exist. 

It is well known that when the cross-sections measured with nuclear targets 
are compared to the free proton results, the two are not identical, but various 
nuclear effects are observed [56l [57l EH] . Although, the QCD factorization 
is not as well-grounded theorem in the case of bound nucleons (see e.g. Ref. 
[59j), the pioneering work [60l[6l] and subsequent analyses like [62l[63], and 
especially the article [IVJ of this thesis, have nevertheless revealed that such 
conjecture holds to a very good precision in describing the world data from 
deep inelastic scattering and Drell-Yan dilepton measurements. In effect, 
only the shapes of the PDFs are modified by the presence of the nuclear 
environment. In other words, although the strong, non-perturbative nuclear 
binding has an effect on the quark-gluon structure of bound nucleons, the 
partons at high appear to largely obey the same QCD dynamics as do 
their free counterparts. 

In what follows, I will describe some general features of global QCD-analyses 
paying special attention to the nuclear PDFs in the light of the publications 
of this thesis. I will keep the discussion here quite condensed, yet logical. 
For much more pedantic description about the free proton global fits with 
comprehensive reference hst, consult e.g. the very profound MSTW paper 
[53] . Also, the lecture notes from the series of Summer Schools ran by 
the CTEQ collaboration are an inexhaustible source of pedagogic up-to- 
date material. Much more technicalities about the nuclear PDF fits can be 
found in the publications of this thesis [Tl | ITTT l HV]. 



7.1 Choice of experimental data 

As mentioned above, the guideline in a global QCD-analysis is to keep it 
really global, i.e. include as many different types of data as possible — 
ideally all. In practice, however, it is necessary to somewhat restrict what 
data is accepted and what is not. 

One obvious restriction is that the factorization framework is only applicable 
when the process is "hard" i.e. the invariant scale inherent for the 
whatever process is large » ^ocd- -^^^ example, in deeply inelastic 
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scattering typical kinematical cuts are 

Q2 > 4GeV^ Ml = iP + qf = + Q^^—^ > UGeV^, 

X 

where the latter condition is to keep away from the resonances. Beyond such 
hmits it may be necessary to account also for the higher-twist Q^^" contri- 
butions (e.g. by parametrizing them as they are usually poorly known), as 
well as target-mass corrections [64} [65]. 

It sometimes happens that independent data sets are not compatible with 
each other. In a case where there are several measurements for the same 
observable and only one data set disagrees with the rest, it may be possible 
to rule out the one measurement as being "wrong" or "not fully understood" , 
and concentrate solely on the others. However, when there are not too many 
independent measurements, it may be necessary to completely abandon that 
type of data for safety i.e. for not being too biased by subjective choices. 
An example of this kind of issue has been the direct photon production 
[66} [671 [68] which is nowadays not included in the latest free-proton PDF 
analyses despite its potential abihty to constrain gluons. However, due to 
complexity of the modern colhder experiments, small mutual inconsistencies 
between independent data sets are more a rule than an exception. This is 
eventually reflected in the PDF uncertainty analysis, necessitating various 
extensions to the strict rules of ideal statistics. 

Typical processes employed in the free-proton PDF analyses include 

• Deeply inelastic scattering related measurements 

• Drell-Yan dilepton production 

• Rapidity distributions in heavy boson (W^ and Z) production 

• Jet measurements 

The sensitivity of these data types for different PDF-components is exten- 
sively documented e.g. in and I will not go to details here. 

Bound protons 

In the case of the nuclear PDF studies the variety, amount and kinematical 
reach of the available data is much smaller, see Fig. I7.1[ For this reason, 
typical nuclear PDF studies have so far adopted the standpoint that the free- 
proton PDFs are taken as fixed, fully known, and only the nuclear modifica- 
tions suggested by the data are inferred. The data conventionally utihzed in 
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10"^ 10"^ 10"^ 10"^ 10"^ 1 

X 



Figure 7.1: Kinematical reach of data employed in typical global analyses 
of the free proton (upper panel) [69j and nuclear (lower panel) PDFs. On 
the upper panel, the kinematical cuts have been already applied, whereas 
on the lower panel the dotted line denote the cut imposed in jlH [nil US 
the DIS data. 

the nuclear PDF studies come not as absolute cross-sections, but as ratios 
of nuclear cross-sections to the free-nucleon ones. In such quantities, things 
hke electroweak corrections and higher-twist contributions can be expected 
to largely cancel, and consequently it is reasonable to push the kinematical 
cuts little further down than is usually done if absolute cross-sections are 
in question. For several years, the only available precision data for nuclear 
PDF analyses were from deeply inelastic scattering and Drell-Yan dilepton 
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production, as shown by the dots in Fig. I7.H but those data alone could 
not constrain the nuclear gluon PDFs very convincingly. Indeed, the usual 
procedure has been to fix a large part of the nuclear modifications for gluons 
by hand invoking reasoning based on physical intuition. To relax the need 
for such assumptions, our latest global analysis in the nuclear sector [Wj 
includes data from inclusive vr'^-production (the gray hne in Fig. 17.10 mea- 
sured recently at BNL-RHIC [70j. These data are not yet very precise, but 
provide a handle to some extent improve the determination of the nuclear 
gluons. 

7.2 The traditional approach 

The conventional procedure of global analyses can be summarized as follows: 
A. 

The PDFs are first parametrized at a cliosen initial scale Qq im- 
posing the sum rules. For absolute free-proton PDFs the functional form 
of the parametrization is typically 

fi{x,Ql)=x-^{l-xfF{x), 

where the function F{x) varies from one analysis to another. Ideally, this 
function should be as flexible as possible, but too much freedom may in- 
duce unphysical, completely artificial features to the PDFs. At the end, 
how much complexity one should build in F{x), depends on the diversity 
and accuracy of the input data. In the nuclear PDF analyses, usually the 
nuclear modification factors Rf{x,QQ) encoding the relative difference to 
the free proton PDFs are parametrized. A typical parametrization indicat- 
ing which x-regions are meant by the commonly used terms: shadowing, 
antishadowing, EMC-effect, and Fermi-motion, is shown in Fig. 17.21 For 
example, in [III Ell HYj we define the bound proton PDFs by 

f,''{x,Q^)^Rf{x,Q'')ff{x,Q^), (7.2) 

for each parton fiavor i. Above f[{x,Q'^) refers to a fixed set of the free- 
proton PDFs, and we consider three different modification factors: Ry{x, Qq) 
for both u and d valence quarks, Rg{x, Ql) for all sea quarks, and Rq{x, Qq) 
for gluons. Also their A-dependence is parametrized but the nuclear cor- 
rections for the Deuteron {A = 2) which are expected to be of order 1..2% 
[62l[71j, are neglected. 
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Figure 7.2: Schematic picture of the nuclear modifications to PDFs. 

In principle there is no real physics reason for assuming flavor-independent 
nuclear modification factors for valence and sea quarks at Qq- Unfortu- 
nately, the presently available experimental data does not allow to deter- 
mine such flavor dependence. However, as was shown in 0, a difference 
between R^^{x, Qq) and R^^{x, Qq) would affect the extraction of the Wein- 
berg weak-mixing angle sin^ from neutrino scattering off an Iron target. 
Concretely, such undertaking was done by the Fermilab NuTeV collabora- 
tion which measured the Paschos-Wolfenstein ratio 172] 



da'^^^ / dxdQ"^ 



da 



u,CC 
Fc 



/dxdQ^ 



^ - sin^ 6'w, 



(7.3) 



where NC and CC refer to the charged (VF^ mediated) and neutral {Z 
mediated) current processes. Surprisingly, the sin^ 6'w found, hed three 
standard deviations away from the world average — an observation known 
as the NuTeV anomaly [73]. Among others sources possible, the simple 
relation above receives a correction from non-isoscalarity of Iron which is of 
the form 



^ - sin^ 6'w 



1 + / 



A 



2Zu^ 



d^ 



A 



+ d6 



(7.4) 



where / is a certain function, and Uy,dy are the average up and down 
valence-quark distributions in a bound nucleon. Thus, a difference between 
i?^^ and would affect the NuTeV analysis. As found in [I], such differ- 
ence can be quite large without really showing up in the global analyses and 
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could very naturally reduce the anomaly by one standard deviation or so. 
There is, however, also correction originating from an asymmetric strange 
sea, J dxx[s{x) — s{x)] / 0, which is not very well constrained either and 
which could explain even the whole NuTeV anomaly alone [7A\ . 

B. 

The absolute PDFs are evolved from the parametrization scale Qq 
to other perturbative scales > Qq by the DGLAP equations. 
On the practical level, an efficient numerical solver for the parton evolution 
is critical for a successful global analysis with a reasonable computing time. 
One, working, solution was discussed in Chapter 4. Another solution would 
be to make a shortcut and use ready-made package hke QCDNUM [75]. 
Under the evolution, also the nuclear modification factors Rf{x, Q^) become 
scale-dependent, and the initial flavour-independence, if it was assumed, 
usually disappears. Significant differences between e.g. R^^ and R^^ do 
not, however, seem to build up if they start from the same initial condition. 

C. 

The cross-sections are computed using the factorization theorem. 
As became exphcitly demonstrated for deeply inelastic scattering (Chapter 
[3]) and Drell-Yan dilepton production (Chapter [4]), the numerical evaluation 
of the NLO cross-sections is sometimes quite demanding, involving mul- 
tiple integrations. Consequently increased computing time is a potential 
"killer" for a global analysis beyond leading order. For example, Monte- 
Carlo integrations for inclusive jet production NLO cross-sections, with var- 
ious experimental cuts implemented, would simply be too slow to be always 
recalculated at every round of parameter iteration. The simplest solution 
is to calculate so-called "K- factors", cjnlo/o'lo, for all data points with an 
educated guess for the NLO PDFs. Then, one can resort to computing only 
leading order cross-section, which, when multiphed by such pre-computed 
factor, serves as an approximation to the NLO one. However, tricks similar 
as used in Chapter [6] to reduce an integration to a matrix multiphcation, 
are also viable [TBI [77]. A special feature of the nuclear PDF studies which 
take the free-proton PDFs fixed, is that sometimes it is possible to perform 
some integrals also beforehand, see Appendix of Ref. [TVj . 

D. 

The computed cross-sections are compared to the experimental 
ones, and the parametrization is varied until the best agreement 
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with the data is reached. What is meant by "best agreement" is, to 
some extent, a matter of convention. An elementary solution is to tune the 
parameters by trial and error to establish a parametrization that simply 
"looks good" [60j when contrasted with the data. However, nowadays such 
procedure is replaced by more algorithmic methods based on minimization 
of a x^-function, defined e.g. as in [Till HV] 



Within each data set labeled by N, Di denotes the experimental data value 
with 5i point-to-point uncertainty, and Tj is the theory prediction corre- 
sponding to a parameter set {a}. The weight factors wn are used to am- 
plify the importance of those data sets whose content is physically relevant, 
but contribution to would otherwize be too small to be noticed by an 
automated minimization. 

In certain cases, an overall relative normahzation uncertainty ^J^""™ is speci- 
fied by the experiment. The normahzation factor /at G [1 - 1 + 5^°™] 
is introduced to account for such cases. For each set of fit parameters {a}, 
its value is solved by minimizing making thus the final not something 
that is put in by hand, but really an output of the analysis. This procedure 
gives a possibility to reproduce the shape (relative magnitude) of the data 
thus capturing the relevant PDF-physics, and not be mislead by problems 
with, sometimes model-dependent, normalization. Apart from the experi- 
mental normalization uncertainty, there may also be unknown contributions 
— hke electroweak corrections — to the absolute computed cross-sections 
having nothing to do with the PDFs. Thus, there is also a theoretical call 
for such tunable normalization factor. 

Although it may sound straightforward, finding the minimum in prac- 
tice is a non-trivial task as the is a highly non-hnear functional of the 
fit parameters. An efficient convergence of a minimization algorithm of- 
ten rehes on the knowledge of the gradient terms and the second-derivative 
matrix at the given location in the space of fit parameters. Due to the 
finite accuracy of the DGLAP-solver and numerical multi-dimensional in- 
tegrations, the becomes actually non-continuous at very small parame- 
ter intervals, which makes the reliable derivative calculations comphcated. 
Consequently, general-purpose packages like MINUIT [78] become easily in- 
sufficient for harsh requirements of global analyses and tailor-made add-ons 
or independent minimizing routines are often needed. 




(7.5) 



N 




(7.6) 
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7.3 Treatment of heavy flavors 



The calculations of the partonic cross-sections a are enormously simplified 
by treating all quarks as massless particles. This is perfectly fine if the 
invariant scales in the cross-sections are much larger than the heavy quark 
masses (only charm and bottom are relevant to present global analyses). 
However, rather than completely forgetting the heavy flavor masses a fol- 
lowing prescription to feed them in is often adopted: 

• The parton distribution for the heavy quark qnix^Q"^) remains zero 
if < m^, but follows the DGLAP evolution when Q'^ > mj^. In 
other words, the number of flavors {uf) is incremented by one every 
time a heavy-quark mass-threshold is crossed. 

• When the factorization scale Q'j in a cross-section is below the heavy- 
quark mass-threshold, the calculation is performed as this quark flavor 
did not exist — even if there is enough center-of-mass energy to phys- 
ically produce it. 

This treatment of heavy flavors is known as Zero-mass variable flavor 
number scheme (ZM-VFNS). It is extremely simple but its shortcomings 
are also rather obvious. In any case, it has been the standard choice in the 
global analyses until only very recently, and it is also the one used in the 
publications of this thesis. 

Another extreme is to work harder with the partonic cross-sections retain- 
ing the full heavy quark mass-dependence and never consider the heavy 
quarks as partons i.e. keep the number of active flavors flxed in the DGLAP 
equations. Such scheme is known as Fixed flavor number scheme (FFNS). 
The problem in this scheme is its hmited domain of apphcability: the par- 
tonic cross-section contain terms ~ log(Q^/m^) which become unstable at 
large Q^, where only resummation of the large logarithms would bring the 
perturbative expansion under control. 

The class of schemes that combine the advantages and avoid the short- 
comings of both ZM-VFNS and FFNS schemes are called General-mass 
variable flavor number schemes (GM-VFNS). In short, these are chains of 
FFNS-type of schemes with specific matching conditions. For a pedagogical 
review, see [79]. These have now become the standard in the most compre- 
hensive QCD analyses and providing a clearly improved fit to the data [80] . 
It is also possible to construct something called Intermediate-mass vari- 
able flavor number schemes [80j, which preserve the simplicity of the ZM 
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scheme, but mimic the full GM scheme by implementing its most relevant 
kinematical effects. 

7.4 Uncertainty analysis 

Finding only the single set of parameters {ao} that optimally fits the ex- 
perimental data in the sense of giving the minimum x^, does not, however, 
alone represent everything that can be learned about the PDFs. Due to 
the experimental uncertainties and fluctuations in the data, the partons 
obtained by steering the fit parameters slightly off from the x^-minimum, 
cannot be right away ruled out as being "completely wrong". Quantifying 
the uncertainties of such origin has become an increasingly important topic 
in the global PDF analyses. For example, the effects of such uncertainties 
in the predictions for the LHC "standard candles", W='= and Z production 
in pp-colhsions, are of the order of few percents [52ll53]. Thus, measuring 
eventually something which is like 10% apart from these predictions would 
be very interesting as such large discrepancy would be difficult to fix by 
simply re-fitting the PDFs. Several methods for performing the PDF un- 
certainty analysis exist. Some common ones and the ideas behind them 
are: 

Lagrange multiplier method |81] 
Defining 

^(A,{a})^x'({a})+AX({a}), 

and minimizing ^'(A, {a}) for several fixed values of A, gives a sequence of 
{x\, Xx)-paLiTS, corresponding to a best that can be achieved with X 
taking a specific value Xx- These pairs comprise a x^-profile as a function 
of X with minimum at X = Xq = X{{ao}). By restricting the allowed 
growth of x^ above the minimum Xoi the achievable range for X becomes 
mapped out. This is a robust procedure, but has a hmited apphcability as 
finding the uncertainty for one single PDF-related quantity X requires an 
abihty to perform several global fits. 

Monte Carlo technique [55J 

The principle is to prepare multitude rephca of the original cross-sections 
cTj by transformation 

ai ai{l + 5iRi) , 

where 5i denotes the experimental uncertainty and Ri is a random number 
drawn from a Gaussian distribution centered at 0. Performing a fit to each 
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of the prepared replica gives a corresponding set of PDFs. The uncertainty 
in any quantity X is then estimated from the spread of the individual pre- 
dictions computed by these sets. Although simple and easy to implement, 
the problem is that in order to get enough statistics, hundreds of separate 
PDF fits are required, but still nothing guarantees that the all relevant 
possibihties are covered ^82] . 



The Hessian method [83J 

This method is superior in its usefulness in various apphcations. As ex- 
plained in detail in [IVj, it rests on expanding the around its minimum 

xl as 

ij * ^ ij 

(7.7) 

which defines the Hessian matrix H. The non-zero off-diagonal elements 
in the Hessian matrix are a sign of correlations between the original fit 
parameters, invalidating the usual error propagation 



for a PDF-dependent quantity X. Therefore, it is useful to diagonahze the 
Hessian matrix, such that 

X'^xl + E'l (7-9) 

i 

where each Zi is a certain hnear combination of the original parameters 
around {a°}. In these variables, the usual form of the error propagation 

is justified. The practicality of the Hessian method resides in constructing 
so-called PDF error sets which are obtained by displacing the fit pa- 
rameters to the negative/positive direction along each Zk separately such 
that x^ grows by a certain amount Ax^. Approximating the derivatives in 
Eq. ()7.10p by finite differences, the error formula can be re-written e.g. as 

(^^f = \E[^(S^)-^(Sk)]\ (7.11) 

A; 

where X{S'^) denotes the value of the quantity X computed with the PDF 
error set St- 
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The methods above require to specify how much the is allowed to gro-wo 
above its minimum value Xo> i-^- what is Ax^. To determine such range is, 
however, far from being a straightforward exercise. On the contrary, it is a 
problematic and much debated issue with no universally agreed procedure. 
As the spirit of the PDF uncertainty analysis is not so much to find statis- 
tically ideal answers, but more to map out the physically relevant range 
for the partons, choosing the statistically ideal one-sigma 68% confidence 
criterion Ax^ = 1 |84] would not make much sense. The basic reason being, 
as mentioned earher, that in a truly global analysis the data sets seldomly 
demostrate a perfect mutual harmony, but small inconsistensies tend to ex- 
ist. Prom the practical point of view, if the absolute minimum is ~ 1000 
(as in the NLO analysis of this thesis [lY]), it would be very naive to rule out 
partons which increase x^ by only one unit. Instead, adopting a prescrip- 
tion that requires each data set to remain roughly within its 90%-confidence 
range [Ml [83l [iVj , we find Ax^ = 50 much better motivated. 

Apart from the experimental uncertainties and their treatment described 
shortly above, there are also other sources of uncertainties. One issue is the 
adopted form of the parametrization which may introduce a bias and make 
the uncertainty bands too narrow in the region which is not covered by the 
data. Such issue can be addressed in a neural network approach, which tries 
to get rid of the fit-function bias [54]. Although providing a very flexible way 
to parametrize the PDPs, it has been criticized as being too flexible: There 
are often physical reasons for picking a fit function of a particular form, and 
too much freedom may lead to unphysical behaviour of PDPs. This is an 
example that falls to the category of "Theoretical uncertainties" [85j among 
other things like unknown higher-order pQCD corrections, resummations 
near the phase space boundaries, higher-twist corrections, and absorptive 
effects in parton evolution. The size of such uncertainties is often impossible 
to quantify in practice. 



^In the Monte Carlo procedure similar role is played by the width of the normal 
distribution. 
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Chapter 8 

Status of global nuclear PDF 
analyses 

Finally, I review the main results from the nuclear PDF analysis of this 
thesis, EPS09 [IV], which is an NLO successor to the pioneering leading- 
order work EKS98 [60j|61j, and articles [II] and [III] of this thesis. 

8.1 Nuclear modifications at LO and NLO 

As mentioned, the goal of these studies is to test the QCD factorization 
and find the process-independent nuclear modifications to the free proton 
PDFs. The results from EPS09NLO are shown in Fig. O for Lead. I 
have plotted the obtained modifications at two scales, at Qq = 1.69 GeV^ 
and at = 100 GeV^, in order to emphasize their scale-dependence. One 
prominent feature that becomes clearly conveyed by this figure is that even 
rather large uncertainty band for the initial small-x gluon modification Rq, 
shrinks along the scale evolution quite a bit. This is a clear prediction of 
the factorization approach. The uncertainty band for the valence quarks 
at large x becomes very small, only ~ 2%, thanks to large amount of deep 
inelastic data. However, as the uncertainty becomes so small, neglecting the 
nuclear effects in Deuteron might not be justified if the PDFs are defined as 
in (17.2p . This is an example of the experimental uncertainties being small 
enough that the theoretical uncertainty may actually dominate. 

The large uncertainty of the small-x gluons is not very surprising. It was 
also investigated in our preceding leading order analysis EPS08 [TTl], where 
we searched for the strongest possible gluon shadowing which was still in 
agreement with the DIS data, but only barely so. Satisfyingly, the result 
found in that analysis is very close to the lower uncertainty hmit which 
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Figure 8.1: The nuclear modifications Ry, Rs, Rg for Lead at an initial 
scale Ql = 1.69 GeV^ and at = 100 GeV^ The thick black lines indicate 
the best-fit results, whereas the dotted green curves denote the individual 
error sets. The shaded bands are the total uncertainty, computed according 
to equation close to Eq. 17.111 



we now find by the Hessian method in the leading-order version of EPS09 
|IV]. This is demonstrated in Fig. 18.21 which displays the low-Q^ nuclear 
modifications for Lead from several leading-order analyses. 



8.2 Data vs. theory 

8.2.1 Deeply inelastic scattering 

The DIS data are the "bread and butter" of all global PDF analyses. In 
Fig. 18.31 1 show a compilation of measured nuclear modifications for deep 
inelastic structure functions with respect to Deuterium 

for various nuclei and compare them with the EPS09NLO parametrization. 
The shaded bands denote the uncertainty derived from the EPS09NLO and 
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Figure 8.2: Comparison of the valence and sea quark, and gluon modifi- 
cations at (5^ = 1.69 GeV^ for Pb nucleus from the leading-order global 
analyses EKS98 ^M\, EKPS P, nDS p], HKN07 [62j, EPS09LO [m], 
and EPS09LO [IV] . 




10"' 10" lo' lo' 10^ lo' 



Figure 8.3: The calculated NLO Rp^{x, Q^) (filled symbols and error bands) 
compared with the NMC 95 [86j and the reanalysed NMC 95 [87] data (open 
symbols). 



— as I would hke to emphasize — their size is very similar to the error bars 
in the experimental data. This a posteriori validates the quite large Ax^ 
which we found as explained in detail in [IVJ. 
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8.2.2 Drell-Yan dilepton production 

The situation is rather similar with the nuclear effects in the Drell-Yan 
dilepton data 



\d(J^/dM'^dxi^2 



rUxi,2, ^ f7p°7^,7;:o7''^ (8-2) 



where is the invariant mass of the lepton pair and xi,2 = \JM'^ j s e^^^. 
1 display the comparison to the E772 and E866 data in Fig. [831 
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Figure 8.4: The computed NLO iJ^yl^) ^^^) (filled squares and error bands) 
as a function of xi (left) and X2 (right) compared with the E866 [89] and 
the E772 [88j data (open squares). 



8.2.3 Inclusive hadron production 

Whereas the free-proton PDF studies can exploit the inclusive jet produc- 
tion data from Fermilab experiments ZO and CDF to access the large- a; 
gluons, there are not yet0 corresponding nuclear data to use for this pur- 
pose. Consequently, the nuclear modifications for the gluons have been 
largely unknown. To do better job in this respect, almost anything which 
is sensitive to the gluons is welcome if it reduces the need for theoretical 
assumptions and fits in with the other data. As we have learned (1111 HY], 

'There is, however, hope that jet measurements performed in nuclear collisions at 
RHIC could be exploited for the nPDF analyses. 
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Figure 8.5: The computed -R^Au (thick black line and blue error band) at 
y = for inclusive pion production compared with the PHENIX [70j (open 
squares) and STAR [90j (open circles) data. The error bars are the statistical 
uncertainties, and the yellow band around PHENIX data indicate the point- 
to-point systematic errors. The data have been multiphed by the optimized 
normalization factor /n = 1.03 for PHENIX and /n = 0.90 for STAR. 



potential candidates are the RHIC measurements for inclusive high-p^ pion 
production, especially for the corresponding nuclear modification 

„^ _ 1 £N^^^/dpTdy niinj,ias ^dP a^^^ / dpT dy 



•-dAu 



(Afeou) d^N^^'/dpTdy d^al^/dprdy 



where the pion transverse momentum is denoted by pt, rapidity by yj?, and 
where (A'cou) denotes the number of binary nucleon-nucleon colhsions. In 
the factorization framework the cross-sections are computed, schematically, 
by 

i,j,k=q,q,g 

where the additional factor D^^.,^ is the fragmentation function for parton 
i to make a pion. Due to the presence of this piece, pion data is not used 
in the free-proton PDF fits as there are uncertainties in disentangling be- 
tween the PDF-originated effects and those inherent to the fragmentation 
functions. However, in the nuclear ratio -RJ^u' ^^'^V details of the frag- 
mentation functions seem to cancel: It is reassuring that the shape of -RJ^u 
is practically independent of the particular fragmentation functions used in 
the calculation — modern sets hke (STJEnillDJ all give about equal results in 
the required kinematical domain [y^ = 0). Having made this observation. 
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it should be quite safe to utilize this type of data in the nPDF analysis. 
The comparison with the PHENIX and STAR data is shown in Fig. [O In 
the computed black curve the downward trend at the large-p-r end is due 
to the presence of an EMC-effect in the large-x gluons, while the turnover 
toward small pt originates from the gluon shadowing. No other effects are 
needed to reproduce the observed spectra, and the fit does not seem to show 
a significant tension between this and the other data. 

As already mentioned, in the EPS08-article [TTl] we studied the chances of 
a very strong gluon shadowing. This study was inspired by the suppres- 
sion observed in the nuclear modification -RdAu for the negatively-charged 
hadron yield in the forward rapidities (r/ = 2.2, 3.2) by the BRAHMS collab- 
oration [91] in d+Au collisions at RHIC. Although a quite deep shadowing 
at small- j; and at low-Q^ was found out to be well within the experimental 
uncertainties, we found a clear tension between the BRAHMS and small-x 
deeply inelastic NMC data, as demonstrated in Fig. 18.61 below. 




Figure 8.6: These three panels stress the tension between the deeply inelastic 
NMC data and the suppression in the forward rapidity hadron production 
-RdAu measured by the BRAHMS collaboration. The three curves demon- 
strate the effect of assigning weights wn = 0, 40, 150 for the BRAHMS data 
in the EPS08 global leading-order analysis: the more weight is given to the 
BRAHMS data, the worse the Q^-slope in the NMC data becomes repro- 
duced. 
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Figure 8.7: Inclusive yields in p+p and d+Au collisions. The experi- 
mental pT-binned data from BRAHMS [91] are shown by open squares with 
statistical and systematic errors added in quadrature. The blue band in- 
dicates the 90% confidence range derived from the CTEQ6.1M and EPS09 
uncertainties. The calculated cross-sections have been averaged over the 
PT-bin width as in the data. 
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Being aware that such issue existed, we decided not to use the BRAHMS 
data in the subsequent EPS09 jiV] analysis. There is, however, a deeper 
reason for not to include the BRAHMS data in EPS09. The Fig. O displays 
the CTEQ6.1M and EPS09 predictions (with fDSS fragmentation functions 
[38] ) compared with the BRAHMS data for the absolute h~ spectra, 

d^NPP minM.s 1 d^gPP _ d^N'^^^ n,in_bias (A^cou) 2^^^^^^^" 

d^pxdy a'^'ff^'^ (Ppxdy ' d^prdy cr^^if^'" d^prdy 

In the r/ = 2.2 panel, the measured p+p and d+Au spectra are both in 
good agreement with the NLO pQCD. However, in the most forward r] = 

3.2 bin there is a systematic and significant discrepancy present between 
the measured and computed pt spectrum for the p+p colhsions, although 
the d+Au spectrum in the same rapidity bin is better reproduced. This 
observation helps to throw some hght on the inconsistency revealed by the 
EPS08-analysis: whether the factorization breaks down already for the p+p 
colhsions and the agreement in d+Au is just by a chance, or there are some 
difficulties with the BRAHMS data in the forwardmost rapidities, where the 
measurements are very challenging. 

8.2.4 Scale-breaking effects 

As the theory goal of the global nuclear PDF studies is largely to test the 
QCD factorization — to find deviations from the DGLAP dynamics — the 
scale-breaking effects deserve special attention. Such effects are cleanly visi- 
ble e.g. in the E886 Drell-Yan data in Fig. 18. 41 where tendency of diminishing 
nuclear effects toward larger invariant mass is observed, although the 
experimental uncertainties are admittedly large. The x-binned DIS data 
versus Q^, shown in Figs. 18.81 and 18. 9| also reveal some general features: At 
small X the Q^-slopes look positive, while toward larger x, the slopes gradu- 
ally die out and become even a bit negative. In both cases, the EPS09NLO 
reproduces such features supporting the vahdity of the factorization. 

8.3 Comparison between the NLO works 

Besides the EPS09NLO, there are two other NLO-level nuclear PDF sets 
available, and Fig. 18. 101 compares the nuclear modifications for Lead as they 
are predicted by these three sets. This comparison is shown at two scales, 
at Ql = 1.69 GeV^ and at = 100 GeV^. The significant discrepancies 
between EPS09NLO and others — that is, curves being outside the blue 
EPS09NLO error bands — are found from the sea quark and gluon sectors. 
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Figure 8.8: The calculated NLO scale evolution (solid black lines and error 
bands) of the ratio F^^/F^ compared with the NMC data |92] for several 
fixed values of x. 
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Figure 8.9: The calculated NLO scale evolution (solid black lines and error 
bands) of the ratio F2 /F2 compared with the NMC data [92] for several 
fixed values of x. 
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Figure 8.10: Comparison of the nuclear modifications for Lead at the scales 
Q2 = 1.69 GeV^ and = lOOGeV^ from the NLO global DGLAP analyses 
HKN07 [62], nDS [63], and EPS09NLO [iVj. 




Figure 8.11: Same as Fig.l8.5[ but also prediction from HKN07 is shown. 
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At low X, the differences rapidly shrink when the scale is increased, but 
at high- a; region notable discrepancies persist. Most of the differences are 
explainable by the assumed behaviours of the fit functions, for which nDS 
and HKN07 allowed less freedom. The use of more rigid form of the fit func- 
tion was due to lack of gluon constraints (they did not employ the RHIC 
pion data which is included in EPS09 [lYj), and giving the gluons more free- 
dom would probably not have improved their fits but instead led to a poorer 
convergence with parameters drifting to their limits. To clarify the conse- 
quences of distinctly different behaviour especially between EPS09NLO and 
HKN07, the Fig. [8lT1 shows the pion i?dAu computed also with HKN07. 
Despite the rather large data uncertainties, the qualitatively contradicting 
sign of the p^-slope in Raku is a promising finding as more data with better 
precision should be soon able to decisively discriminate between different 
scenarios. 

8.4 Conclusions and future prospects 

At the end, the conclusion that can be drawn from the global QCD analyses 
for nuclear PDFs performed so far is that the QCD factorization conjecture 
seems to work very well in the explored kinematical region and for the 
presently analysed data types. In order to find evidence of discrepancies — 
especially nucleus-enhanced power corrections in the parton evolution [93] 
or in the cross-sections [94] — the scope of the global analysis should be 
therefore extended. In addition to enlarging the kinematical coverage and 
the variety of included processes, also an extended theoretical treatment 
of the parton dynamics (by including power-suppressed terms) should be 
of importance. The next steps at the horizon toward reaching smaller x 
and higher Q^, could possibly be realized by running the LHC also in the 
p+Pb mode. In farther future, electron-ion colhders hke the planned eRHIC 
[95] or LHeC [96], would also penetrate deeper in the (x, Q^)-plane. About 
new processes to include, direct photon data from RHIC d+Au and Au-Au 
colhsions should be available shortly and there are also data from deeply 
inelastic neutrino-Iron scattering from the Fermilab NuTeV collaboration 
[97] available. The neutrino data has recently been claimed [98] to point 
toward different nuclear effects than what are obtained e.g. in EPS09. If 
true, it would indicate that there are process-dependent effects present, 
casting doubts on the factorization. This issue certainly deserves a further 
clarification. 
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It is a sort of a peculiar situation that almost all global free-proton fits 
include data with nuclear (deuterium or heavier) target^ — all groups with 
their own favorite way to "correct" for the nuclear effects. Thus, there is 
a danger of circularity as the basehne PDFs which should be free of any 
nuclear effects, do actually somewhat depend on how the nuclear effects 
were corrected for — the free-proton and nuclear PDF analyses are not 
entirely independent. This observation opens the road for a future work 
which will combine the free and bound proton PDF analyses into a single, 
"master" global analysis. 



^For example, the NuTeV data is one of the main sources for constraining the strange 
sea asymmetry. 
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